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PREFACE 


This  report  discusses  the  motion  of  gases  and  weapon  debris  and  the 
temperature  distribution  in  the  immediate  vicinity  of  the  rising  late- 
time  fireball  of  a  nuclear  detonation  in  the  atmosphere.  "Late-time" 
loosely  denotes  that  period  after  the  fireball  has  adjusted  to  the  pres¬ 
sure  of  the  ambient  atmosphere.  This  subject  will  be  of  interest  to 
those  agencies  which  seek  methods  for  determining  the  influence  of 
the  tireball  on  electromagnetic  wave  propagation. 

The  authc  '  wish  to  acknowledge  the  assistance  of  R.  W.  Hendrick, 

T.  J.  Barrett,  G.  Schindler,  H.  Kramer,  E.  Tschupp,  and  G.  Rohringer 
of  General  Electric,  TEMPO,  who  have  entered  into  a  number  of  dis¬ 
cussions  on  the  subject  of  this  report. 

The  co-author,  Dr.  Charles  Dolph,  is  Professor  of  Mathematics  at 
the  University  of  Michigan  and  consultant  to  the  General  Electric 
Company . 
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SUMMARY 


KM  cSTMP-/ 


This  report  describes  a  teclinique  for  the  solution  of  a  class  of  prob¬ 
lems  in  compressible  turbulent  flow  and  demonstrates  its  application 
to  a  problem  involving  the  rise  and  the  entrainment  of  ambient  air  by 
the  late -time  fireball  due  to  a  nuclear  detonation  in  the  atmosphere. 

The  method  is  based  on  the  work  of  J.  E.  Moyal'  which  shows  that  the 
components  of  Fourier  spectra  of  velocity  parallel  to  the  wave  vector 
are  associated  wivh  compression  waves  (irrotational  flow),  while  those 
components  transverse  to  the  wave  vector  are  associated  wTth  eddy 
turbulence.  If  solutions  are  sought  for  systems  in  pressure  equilibrium 
with  an  undisturbed  atmosphere,  the  resulting  formalism  takes  on  con¬ 
siderable  simplicity  due  to  uncoupling  of  the  equations  of  motion 

Assumptions  related  to  th  •  evaluation  of  certain  convolution  integrals 
which  appear  in  the  theory  are  examined  and  shown  to  involve  the  larg¬ 
est  eddies  In  the  role  of  the  principal  convective  agent.  This  assump¬ 
tion  converts  the  nonlinear  Navier-Stokes  equation  to  liner  r  form  with 
solutions  conveniently  described  in  a  moving  frame  of  reference.  The 
heat  transfer  and  continuity  equations  are  also  simplified  by  the  same 
linearization  process.  Simultaneous  solutions  of  the  equations  in  wave- 
vector  space  and  subsequent  inversion  to  physical  space  v.dth  certain 
initial  conditions  provides  solutions  which  have  considerable  similarity 
to  the  actual  flow  situation  fround  in  weapons  tests  such  as  vortex  ring 
(torus)  formation.  Qualitative  and  quantitative  features  of  the  predicted 
motion,  for  exp.mple ,  velocity  distribution,  rise  rate,  and  temperature, 
are  discussed. 

A  section  on  conclusions  and  recommendations  assembles  the  main 
ideas  and  results  which  have  come  from  the  study  and  suggests  the  di¬ 
rection  of  future  effort  both  with  respect  to  the  fireball  problem  and 
with  respect  to  the  general  theory  of  turbulent  flow. 

A  set  of  appendices  is  supplied  which  provides  information  on  the  lit¬ 
erature  in  the  field  of  turbulent  flow,  a  method  of  solving  the  linearized 
equations  of  Moyal  for  large  Reynolds  number,  and  provides  inf  irma- 
tion  on  certain  integrals  used  in  the  text. 
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SECTION  1 

INTRODUCTION 


A  task  on  "Investigatibn  of  the  Effects  of  Nuclear  Detonations  on  Elec¬ 
tromagnetic  Propagation,"  Contract  AF  19(628) -4210,  has  been  the 
study  of  turbulent  mixing  processes.  The  objective  of  this  s'^udy  has 
been  to  obtain  information  on  the  mixing  processes  in  the  fir":hail 
which  might  determine  the  distribution  of  debris,  temperature,  and 
velocity. 

The  motion  and  temperature  distribution  of  the  late-time  fireball  has 
considerable  importance  from  the  point  of  view  of  the  weapons  system 
engineer.  This  is  related  to  at  least  two  separate  areas  of  penetration 
aid  and  ABM  system  design:  the  absorption  of  EM  radiation  (such  as 
ECM  and  radar  waves)  which  is  dependent  on  temperature  distribution 
and  debris  motion;  and  clutter  which  is  due  to  irregularities  of  the  prop¬ 
erties  of  the  fluid.  Both  of  these  effects  strongly  influence  electronic 
system  design.  Because  of  this  importance  a  number  of  investigations 
of  the  phenomenon  associated  with  the  rise  and  entrainment  of  the  fire¬ 
ball  have  been  conducted.’!*  Test  information  h?s  been  correlated  and 
empirical  expressions  developed  for  average  temperature,  rise  veloc- 
ity.  and  fireball  radius.  These  expressions  are  useful  kinds  of  infor¬ 
mation  for  the  system  designer  since  over  certain  ranges  of  yield  and 
altitude  of  detonation  they  orovide  data  for  ^he  gross  features  of  the 
fireball  such  as  an  average  terr.perature .  The  variations  of  temper¬ 
ature  and  the  detailed  m.otions  in  the  vicinity  of  the  rising  fireball  have 
not  been  forthcoming  up  to  now  largely  because  the  problem  involves 
the  complicated  concepts  of  gravitational  convection  in  a  compressible, 
viscous,  heat  conducting  fluio. 

This  report  attempts  to  provide  some  further  information  on  the  fire¬ 
ball  rise  problem  by  using  the  most  general  methods  of  fluid  mech.anics. 
The  problem  is  first  posed  in  as  complete  a  form  as  possible.  Through¬ 
out  the  analysis  a  series  of  assumptions  is  invoked  to  continuously  sim¬ 
plify  the  problem  in  a  way  which  preserves  the  main  mechanisnis  of  the 
dynamics  and  still  allows  solutions  to  be  obtained  which  display  some 
important  features  of  the  problem. 

’i'See  References  2,  3,  4  and  5. 
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Unfortunately  the  problem  in  its  fullest  extent  must  be  posed  within 
the  framew’ork  of  compressible  turbulent  flow'.  There  has  been  only 
a  small  amount  of  fundamental  work*  (see  Appendix  B)  in  this  field 
although  applications  involving  mixing  length  or  eddy  viscosity  as  ' 
sumptions  in  wakes,  jets,  and  supersonic  boundary  layers  are  highly 
developed.  The  work  on  turbulent  gravitational  convection  has  involved 
assumptions  such  as  similarity^®  or  restrictions  on  the  x’ole  of  density 
differences , 

Although  much  important  insight  into  the  problem  has  been  made  avail¬ 
able,  by  this  means  it  has  been  difficult  to  find  full  solutions  which  yield 
the  kinds  of  motion  w'hich  experience  shows  the  fireball  to  undergo  such 
as  the  formation  of  a  moving  vortex  ring.  Because  of  this  it  appeared 
that  a  different  class  of  assumptions  would  lead  to  more  realistic 
motions . 

O.  M.  Phillips  has  shown  that  if  the  net  linear  momentum  of  the  fluid 
is  non-zero  in  a  viscuous  fluid  in  the  final  stages  of  decay  of  turbulence, 
the  motion  develops  into  a  kind  of  vortex  ring  similar  to  the  solutions 
found  in  this  report.  The  basic  difference  is  that  the  inclusion  of  part 
of  the  non-linear  term  as  carried  out  here  provides  a  translation  of  the 
ring  through  the  fluid  which  depends  on  the  initial  conditions.  Phillips 
demonstrates  that  the  low  density  in  the  core  of  a  heated  vortex  ring 
may  be  expected  to  keep  the  material  of  the  ring  coherent,  presumably 
because  of  the  effect  of  centrifugal  forces  operating  on  the  density  dif¬ 
ferences.  The  theory  of  vortex  rings  in  an  inviscid  fluid  due  to  Lamb 
and  to  Hil)  are  reviewed  by  Turner.  Some  excellent  photographs 

of  rings  are  contained  in  Reference  14. 

In  an  important  paper  Moyal^  demonstrates  the  inethodology  by  which 
two  types  of  velocity  spectra  may  be  segregated.  The  resulting  equa¬ 
tions  are  tractable  if  the  a",sumption  is  made  that  the  largest  eddies 
in  the  flow  (the  only  ones  which  couple  wHh  the  constant  gravitational 
field)  dominate  the  non-linear  terms.  This  last  assumption  appears 
reasonable  when  the  equations  of  the  velocity  spectra  are  examined. 

The  convolution  integrals  which  arise  from  the  non-linear  terms  after 
Fourier  transformation  are  approximated  by  the  value  of  the  integrand 
which  corresponds  to  non-linear  coupling  between  the  largest  eddies 
and  the  eddies  of  all  other  sizes.  It  is  shown  that  the  resulting  term 
is  responsible  for  translation  of  the  space  coordinates  in  a  manner 
which  depends  or  che  initial  conditions.  The  translation  effect  allows 

*See  References  6,  7,  8  and  9- 
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a  convenient  description  of  the  motion  and  temperature  distribution  in 
a  moving  frame  of  reference.  For  the  case  of  an  initial  velocity  in  the 
upward  direction  only,  the  translational  velocity  may  be  identified  with 
the  rise  velocity  of  the  fireball. 
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SECTION  2 

THE  EQUATIONS  OF  MOTION,  HEAT  TR  NSFER, 
AND  STATE  IN  PHYSICAL  SPACE 


The  problem  is  specified  in  terms  of  the  conservation  of  mass,  mom¬ 
entum,  energy,  and  the  equation  of  state  using  the  dependent  variables 
s  ,  TTp  ,  and  9  where  s--fn{p/p^),  n  =-fn{p/pJ,  B=ln{T/T^)  and  the 
velocity,  v  =  {vj  ,  ,  v:^)  of  the  fluia.  The  subscript  zero  corresponds 

to  any  convenient  time -independent  reference  condition  and  p,  p,  and 
T  are  the  density,  pressure  and  temperature  at  any  point  x'=  {xj  ,  ,  X3  ) 

and  any  time  t^ .  We  specify  a  body  force  per  unit  mass  of  g  in  the 
negative  x'^  direction. 
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The  summation  convention  is  used  in  these  eouations  and  throughout 
this  report.  In  the  order  given,  the  equations  represent  the  conser¬ 
vation  of  mass,  momentum,  and  energy  and  the  equation  of  state. 
The  quantities  p,  x  are  the  kinematic  viscosity  and  thermal  diffus- 
ivity  of  the  fluid,  respecti  ely.  The  quantities  a',  y,  c  are  the 
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velocity  of  sound,  ratio  of  specific  heats,  and  the  specific  heat  at  con¬ 
stant  pressure  of  the  fluid,  respectively.  The  quantity  6^3  is  equal  to 
unity  if  i  =  3  and  zero  otherwise.  In  performing  the  following  analysis 
we  make  the  assumption  that  the  quantities  V,  x,  a',  y,  and  Cp  are 
constant. 


In  many  problems  involving  the  late -time  fireball,  with  low  Mach  num¬ 
ber  rise  rate,  the  term  in  braces  on  the  right-hand  side  of  Equation  3, 
representing  the  effect  of  viscous  dissipation  on  the  temperature,  is 
small  compared  to  the  other  terms  and  may  be  disregarded.*  This 
does  not  imply  that  viscous  dissipation  is  neglected  insofar  as  it  is  the 
ultimate  energy  sink  in  the  smallest  eddies.  It  is  simply  not  consid¬ 
ered  to  influence  the  temperature  measurably. 

In  addition,  the  equation  of  state  represents  a  perfect  gas  and  we,  there¬ 
fore,  assume  that  the  effect  of  chemical  reactions  within  the  fluid  plays 
a  negligible  role.  This  contributes  to  inaccuracy  of  the  solutions  when 
the  temperature  is  above,  say,  2000°K  at  most  altitudes. 

It  is  implicitly  assumed  in  writing  the  above  equations  that  the  flow  be¬ 
haves  according  to  continuum  mechanics. 

We  have  also  assumed  that  the  fireball  is  small  compared  to  distances 
over  which  the  earth's  gravitational  field  chav^es  by  a  significant  amount. 

The  late-time  fireball  is  characterized  by  a  relatively  uniform  pressure 
at  a  given  altitude.  This  concept  of  "pressure  equilibrium,  "  though  it 
is  probably  not  completely  valid,  will  be  employed  to  simplify  the  equa¬ 
tions.  This  assumption  requires  that  the  pressure  everywhere  in  the 
fireball  is  equal  to  that  which  occurs  in  an  undisturbed  atmosphere. 

This  requires  that 
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where  H'  is  the  scale  height  of  pressure  in  the  atmosphere  and  is 
given  by 


and 


If  Equations  5  and  6  are  employed  in  the  momentum  and  energy  equa¬ 
tion  these  become 
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The  other  equations  remain  unchanged.  Equations  1,  4,  8,  and  9  are 
toe  basic  set  which  will  be  discussed  in  the  following  sections. 
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SECTION  3 


THE  NON-DIMENSIONAL  EQUATIONS 
IN  WAVE -VECTOR  SPACE 


Following  the  work  oi  Moyal  the  set  of  Equations  1,  4,  8  and  9  is  trans¬ 
formed  into  wave-vector  space  in  accordance  with  the  following  relations 


dT<:  p  -  i  k  •  X 

dZ.  = - 7  1  V.  e 

(277) 


dx 


+ik • x' 


dZ. 


(10) 


The  quantity  dZ-  is  the  Fourier  spectra  of  the  velocity  in  the  i  di¬ 
rection.  The  integral  on  the  right  is  shown  in  Stieltjes  form  to  allow 
for  the  possibility  that  dZ^  may  be  a  "oathological"  function  (or  im¬ 
proper  function  in  the  sense  of  the  Dirac  delta  function). 


Similarly,  we  define  formally  the  quantities  dS  ,  dP  ,  and  dQ  as 
follows: 


dS  = 


dk 


(2T7) 


-  i  k  •  x'  P  +  i  k  •  x'  ^ 

s  e  dx  ;  s  =  !  e  dS 
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ree-“'’*'  di'  ;  8=  ’‘'dQ  . 


The  integrals  on  the  left  are  taken  over  all  of  physical  space  while 
those  on  the  right  are  taken  over  all  of  wave- vector  space. 

The  magnitude  of  the  vector  k  =  (k^  ,  k^  .  k3)  can  be  loosely  consid¬ 
ered  as  an  inverse  measure  of  eddy  size  in  the  sense  that  277/eddy 
lengtl  i  the  i  direction  «  k^  .  The  transformation  into  wave-vector 
space  thus  brings  together,  into  a  small  volume  in  this  space,  all 
eddies  in  the  flow  of  a  given  size  and  shape  independent  of  where  they 
may  occur  in  physical  space.  The  largest  eddies  in  the  flow  are 
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transformed  so  that  they  reside  near  the  origin  in  wave-vector  space, 
w'hile  the  smallest  eddies  occupy  regions  far  from  the  0'*igin  in  that 
space.  Thus,  within  the  concept  of  flow  of  kinetic  energy  in  w\ave- 
vector  space,  energy  is  considered  to  start  with  the  largest  eddies 
near  the  origin  and  flows  outward  tow'ard  the  region  occupied  hy  the 
smallest  eddies  where,  by  the  action  of  viscosity,  it  is  eventually  dis¬ 
sipated  as  heat.  In  this  manner  energy  is  removed  from  the  gross 
motion  of  the  fluid  as  the  eddies  are  reduced  in  size  and  reappear 
more  remote  from  the  origin.  This  process  cannot  continue  indef¬ 
initely  without  supplying  energy  to  the  flow,  mainly  to  the  larger  ed¬ 
dies.  The  role  of  the  body  force  per  unit  mass  in  Equation  Z  is  to 
supply  this  energy  from  the  gravitational  field.  This  concept  will  be 
used  later  in  an  attempt  to  reduce  to  manageable  form,  the  convolution 
integrals  which  will  soon  appear. 

The  transformation  is  straightforward  and  the  non-dimensional  e  jua- 
tions  in  wave -vector  space  are 


Re  k.  dZ.  -  i  Re  I  k'  dZ  .  (k' -  k)  dS(k') 
J  J  J  J 


(11) 


^dZ.  +  k^  dZ  +  -  k  k  dZ  -  -  i  Re  "k'.  dZ.{k'-k)  dZ.{k')  .  (IZ) 
St  1  1  3  i  J  J  -  J  J  1 


^  dQ  +  -  k"  dQ  =  -  iRe  i'k'.  dZ.(k'-  k)  dQ(k') 
dt  V  •-  J  J 


—  i  Re 
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dP  =  dS  +  dQ 


(13) 

(14) 


Adding  Equations  11  and  13  we  have  a  useful  auxiliary  equation 


■y/  2.  R  C 

y  ~  k  dQ  -  i  yRe  k.  dZ.  ^  dZ, 

t'  J  J  H  3 


(15) 


In  these  equations  a  typical  length  in  the  flow,  r^  ,  is  taken  as  the 
unit  of  length;  a  typical  velocity,  v^  ,  becomes  the  unit  of  velocity. 
The  Reynolds  number,  Re  ,  is  given  by  r^v^/p  .  The  unit  of  time 
becomes  fg/t'  •  Therefore,  all  quantities  in  the  set  of  equations  im¬ 
mediately  above  are  non-dimensional,  including  k  which  is  measured 
in  units  of  1  / r 
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SECTION  4 

DECOMPOSITION  OF  THE  VELOCITY 


Moyal  has  shown  that  the  component  o£  the  vector  dZ  in  wave-vector 
space  which  is  parallel  to  the  vector  k  is  associated  with  compression 
w'aves  and  random  noise  (ir rotational  flow).  The  components  of  dZ 
which  are  transverse  to  k  are  associated  with  the  eddy  turbulence 
(solenoidal  flow).  In  order  to  segregate  these  two  kinds  of  motion, 
following  Moyal,  we  define  dZj^  as  the  component  of  dZ  along  k  and 

dzj'”^  as  the  ith  component  of  dZ^*-)  a  vector  transverse  to  k  .  This 
may  be  written 

dZ.  =  dz:’^'  +  —  dZ,  .  (16) 

1  1  K  k 

Note  that  k-  dZ-  =  k  dZ.  and  kj  dZj*"^  ■  0  . 

If  we  define  the  quantity  dCj(X)  by 

dC.(X)  =  -iRe  r  k'  dZ  (k' -  k)  dX.(k')  . 

^  k'  J  j 

then  we  have 


dC.(X:  ^  dC:^^  (X)  +  dC  (X)  , 

X  1  K 


where  dX  is  any  one  of  the  spectra. 

If  Equation  16  is  substituted  into  Equation  12  and  both  sides  are  multi¬ 
plied  by  kj^/k  the  result  is 


^  dZ^  .  I  k'  dZ^  =  dC,^(Z) 


If  Equation  12  is  multiplied  by  (6..-k.k./k  )  the  result  is 

ij  1  J 
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3t 


=  dci^^  (Z) 
J 


(18) 


The  quantity  ^ij-k^kj/k^  is  called  the  projection  operator  and  has  the 
property  of  converting  a  vector  in  k-space  into  components  normal  to 
k  .  That  is,  if  Aj  is  a  vector,  the  vector  Aj  is  normal  to  k  when 

^  =  A.  (6  .  -k.k./k^  )  . 

1  J  \  ij  1  J  / 

We  now  have  two  equations  which  result  from  the  original  momentum 
equation.  One  of  these  equations  governs  the  behavior  of  the  spectra 
of  the  irrotational  velocity  (Equation  17);  the  other  governs  the  be¬ 
havior  of  the  spectra  of  eddy  turbulence  (Equation  18). 
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SECTION  5 
LINuARlZATION 


The  non-linear  effects  represented  by  the  convolution  integrals  will 
be  partially  taken  into  account  by  assuming  that  the  largest  eddies  in 
the  flow  (those  for  which  k  =  0)  are  the  principal  agents  in  the  convec¬ 
tive  process.  It  is  possible  to  decompose  the  convolution  integrals 
into  the  following  form; 


Fk'  dZ  (k'-  k)  dZ.(k') 
-  e  e  j 


f 

kVk 


k'  dZ  (k'-k)  dZ.(k') 
e  e  j 


k  k 


1 

k 


and 


+  k  dZ  (0)  dz!^^  (k) 
e  e  1 


k'  dZ  (k'-k)  dZ.(k')  k  /k  =  f  k'  dZ  (k' -  k)  dZ  .(k')  k./k 
e  =  J  J  =  e  J  J 


+  k  dZ  (0)  dZ,  (k) 
e  e  k 


where  we  have  used  the  relations 

k.  dZ.  =  kdZ,  . 

1  1  k 


dz!^^k) 

1 


(6. .  -  k.k. /k^  I  dZ  .(k)  and 
\ij  1  J  /  y 


The  last  term  on  the  right  in  these  equations  represents  the  non-linear 
interaction  of  the  largest  eddies  with  the  eddies  of  all  other  sizes  and 
is  the  part  of  the  integral  which  corresponds  to  the  spherical  shell  in 
wave-vector  space  for  which  k'"k  .  The  non-linear  effects  associated 
with  the  first  term  on  the  right-hand  side  are  disregarded  for  the  pres¬ 
ent  study.  Only  the  largest  eddies  are  coupled  to  the  buoyant  force. 
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The  equai'ions  for  the  spectra  of  velocity  become 


3  (t)  2  (t)  (t) 

T-  dz:  '  +  k  dz:  =  -  i  Re  k  dZ  (0)  dZ)  '  (k) 
ot  1  1  e  e  1 


^dZ,  +  I  k"  dZ.  =  -iRek  dZ  (0)  dZ,  (k) 
ot  k  3  k  e  e  k 


The  time  dependence  of  the  quantity  kg  dZg(O)  is  found  by  letting 
k=0  in  the  above  equations 


^  dz!^^  (0)  =  0 

St  1  ' 


(19) 

(20) 


AdZ^(0)=0  . 

The  equations  may  be  integrated  to  give 

dZ^*^^  (0,  t)  =  dz!^^  (0,  0) 

1  1 

dZ^(0,  t)  =  dZ^(0,  0)  . 


The  resulting  inhomogeneous  equations  for  the  spectra  of  the  solenoidal 
velocity  are  linear*  where  the  relation  kgdZj^=kdZj^  has  been  used  again 


+  i  Re  k  dZj^(0 


.0) 


dz!’^^  =  0 
1 


The  equation  for  the  spectra  of  ^  .e  irrotational  velocity  is  identical 
except  that  the  term  k^  is  replaced  by  4k^/3  . 


The  solutions  are 


dz!'^^  (k,  t)  =  dz|*^^  (k,  0)  e  exp  i  Re  k  dZj^(0,  0)  t 


‘] 


(21) 


*If  the  bouyant  force  term  is  not  equated  to  the  pressure  gradient  and 
carried  along  to  this  point  the  resulting  linear  equations  have  time- 
dependent  coefficients. 
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SECTION  5 


'  3  ^  r  1 

dZ  (k,  t)  =  dZ,  (k,  0)  e  exp  -  i  Re  k  dZ  (0,  0)  t 

k  k  L  K  J 


(22) 


Note  that  dz!^^  (k.  0)  =  (6. .  -  k  k . /k^  jdZ  ,(k,  0) 
1  \  3  J  /  J 


k  dZ,  (k,  0)  =  k.  dZ.(k.  0)  . 

k  J  J 


The  total  velocity  spectra  are  from  Equation  16: 


m 


dZ. (k,  t) 


dz|^^  (k.  t)  +  (k./k)  dZj^(k,  t) 

=  dZ.(k,  0)  exp  i  Re  k^  dZ^(0,  0)  tj 


)..k^’  -k.k.\  - 
1'  1  J/ 


2 

-  k  t 

—  +  k.k. 
1  J 


■  I’'  ‘ 


(23) 


The  first  factor  involv'ing  the  exponential  translates  the  resulting  in¬ 
version  in  physical  space  so  that  each  xj  is  replaced  by  x.  -RedZj^ 
(0,  0)t  .  The  coordinate  system  in  which  the  motion  is  most  conveni¬ 
ently  described  is  one  which  is  moving  with  non-dimensional  velocity 
RedZj(0,0)  in  each  of  the  Xj  directions.  The  factors  6— k^'-k.k. 
and  k^k  become  differential  operators  on  the  inverse  of  the  functions 

**2  2 
-  k  t  ?  -  4/3  k  t  2 

dZj(k,  0)  e  /k^  and  dZj(k,  0)  e  /k  .  As  will  be  shown 

later,  for  certain  initial  conditions,  the  differential  operators  cause 
vortex  and  torus  (or  "smoke  ring")  formation. 

If  Equations  13  and  15  are  combined,  the  energy  equation  becomes* 

—  dQ  +  y  -  k^  dQ  =  -  i  Re  k.  dZ.(0,  0)  dQ  -  i(y-  1)  Re  k  dZ. 

u  J  J  J  J 


(24) 


^Attempts  to  solve  the  continuity  equation  for  density  yielded  integrals 
which  could  not  be  evaluated. 
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witu  solution 


dQ(k,  t)  =  exp  (-  i  Vt 


vt->  ik’tj  j 


nt 


dQ(k,  0) i(y-  i)  Re  :  k,dZ.(k.t') 
(  o  J  J 


X  2 

exp  -  i  Vt'  +  y  —  k^*  t  dt 

L  ^ 

where  we  have  defined  V  by 


V  -  Re  k^  dZ^(0,  0)  . 

From  Equation  23  with  the  translating  term  as  a  factor  we  have 


4  ,  - 

-  -  k  t 


k.  dZ.{k.  t) 
J  *1 


So  that 


=  k.  dZ.(k.  0)  e 
«J  *1 

exp  (  -  i 

X  ,  2 
y  —  k  t 

xp  (  -  i  Vt)  ■ 

dQ(k.O)  e 

V 

/  4  , 

2\  / 

X  ,2  \ 

exp  -  k 

•1  -exp(- 

y-  k  ,) 

J  J 


(>  5  -  5) 


(25) 


To  insure  compatible  initial  conditions  of  velocity  ana  temperature 
we  use  Equation  1  5  in  the  form 


dQ(k,  0)  =  -  i  y  kA  dZ.(k,  0) /k'^ 

y  21  V  J 

u 


(26) 


substituting  into  Equation  25  we  have  for  the  spectrum  of  the  temp¬ 
erature  function 
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rpV'M  J 


dQ(k,  t)  =  Re  dZ^(k,  0)  exp  (  -  i  Vt) 


i{?-  l)k. 

^  Hi) 


X  ,  I 

y  -  k  t 

V 


Pr 

yH 


-  i  Pr  k 


J 


4  2 

i(y-  l)k.  “  3  ‘ 

_  J  f _ 


H-l) 


where  Pr  =  y/x  .  the  Prandtl  number. 


(27) 


Again,  the  term  exp(-iVt)  indicates  that  a  moving  coordinate  sys¬ 
tem  is  the  convenient  frame  of  reference  for  a  description  of  the 
temperature. 


This  will  be  inverted  in  the  next  section  for  a  specific  initial  velocity 
distribution. 


The  density  function  may  be  found  directly  from  the  temperature  func¬ 
tion  and  the  equation  of  state. 
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SECTION  6 


QUANTITATIVE  FEATURES  OF  THE  SOLUTIONS 
FOR  THE  FIREBALL  PROBLEM 


SELECTION  OF  INITIAL  CONDITIONS 

•Shortly  after  the  fireball  has  come  to  pressure  equilibrium  it  is  moving 
upward  due  to  ballistic  and  buoyant  forces.  The  velocity  is  likely  to  be 
predominately  in  the  upward  direction  and  to  be  largest  at  the  center 
and  reducing  in  magnitude  toward  the  edges  of  the  visible  fireball.  We 
assume  that  the  initial  motion  can  be  represented  by  a  single  compo¬ 
nent  of  velocity  in  the  z(or  x^)  direction.  The  upward  initial  velocity 
is  assumed  to  follow'  a  Guassian  distribution  so  that  the  initial  conditions 
are 


v^(x,  0) 


30 


r^  /  4b 


Vj(x,  0)  =  ^^{x,  0)  =  0 


(28) 


w’here  v^q  is  the  upward  velocity  at  the  center  and  b  is  a  parameter 
governing  the  rate  at  which  the  velocity  reduces  toward  the  edge  of  the 
fireball.  The  origin  of  the  stationary  coordinate  system  now  requires 
definition  and  is  taken  at  the  center  of  the  fireball  at  t  =  0.  The  quan¬ 
tities  r  and  b  are  non-dimensional. When  r  ”2vb  the  velocity  has 
dropped  by  a  factor  of  e  below  its  value  at  the  center  of  the  fireball. 
Transformation  of  the  initial  conditions  gives 


dZ3(k.  0) 

dZ3(0,  0) 


b  j 

3/2  -k^b 

(29) 

^30  * 

1  e 

i  b  \ 

3/2 

^^30! 

lil 

• 

(30) 

THE  VELOCITY  FIELD 

The  inversion  of  Equation  23  with  the  initial  conditions  specified  above 
in  a  frame  of  reference  moving  upward  with  velocity  Rev^Qlb/ff)^/^  is 
(see  Appendix  C) 

’^Normally  we  will  choose  b  =  1/4  so  that  the  velocity  has  reduced  by  a 
factor  e  when  r=  1. 
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where  S(u)  =  uerf^u—  erfu  . 

The  functions  T(u)  and  S(u)  are  plotted  versus  u  in  Figure  1. 
Note  that  as  t-0  ,  ^^-*72  =  T  and  we  have 

Vj(x,0)  =  v^Cx.O)  =  0  ,  t  -  0 

v^{x,0)  =  -2/i  erf''  ,  t 0  ■ 

T  r 

The  derivatives  of  the  error  function  are 


(34) 


./  2  — u 

erf  u  =  - e 

/? 


t*  • 

erf  u 
Yielding 


JH 


2 


— u 

e 


V3{x,0) 


30 


-r^74b 

e 


The  initial  conditions  have  been  recovered  at  t=  0. 


The  horizontal  radial  velocity  Vp  given  by  Equation  33  is  seen  to 
change  signs  as  we  move  from  above  to  below  the  origin  of  the  moving 
coordinate  system.  It  is  zero  on  the  vertical  axis  and  on  the  horizontal 
plane  through  the  moving  origin.  It  rapidly  approaches  zero  as  r-*'®  ■ 

The  general  character  of  the  horizontal  motion  is  as  shown  in  Figure  2, 
where  the  velocity  is  shown  outward  from  the  vertical  axis  above  the 
horizontal  plane  containing  the  moving  origin  and  tow'ard  the  axiS  below 
this  plane. 
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The  vertical  velocity  has  a  more  complex  dependence  on  time  and  the 
space  coordinates.  In  the  moving  frame  of  reference  it  is  symmetrical 
about  the  plane  =0,  positive  near  the  vertical  axis,  and  negative  ^ar 
from  the  axis  for  small  values  of  time.  As  time  increases  it  reduces 
in  value  near  the  vertical  axis  and  increases  at  values  of  p  which  are 
increasingly  large.  Its  behavior  is  as  shown  in  Figure  3. 


When  combined  with  the  horizontal  radial  velocity  Vp  a  configuration 
in  the  moving  frame  shown  in  Figure  4  results.  This  displays  the  fea¬ 
tures  often  observed  in  nuclear  weapons  tests  (and  chemical  explosions) 
where  a  toroid  or  vortex  ring  is  formed. 


MOVING 


EARLY  TIME  UTE  TIME 

Figure  3.  Configuration  of  the  upward  velocity  in  the  moving  frame  of  reference 
in  the  plane  X|  =  0. 


TIME  - ^ 

Figure  4,  Velocity  configuration  at  three  different  times  showing  vortex  ring 
formation  in  a  frame  of  reference  moving  upward  v/ith  uniform 
velocity. 


20 


SECTION  6 


TEMPERATURE  DISTRIBUTION 


If  we  continue  to  employ  an  initial  upward  velocity  which  has  a  Gaussian 
distribution  in  space  about  the  center  of  the  early  fireball,  the  tempera¬ 
ture  spectra  is  given  by 


3/2 

dQ(k,t)  =  Re  v^^  (b/ w)  expL-iVtJ 


/ 

I  El 

\yH  ~ 


i(y-l)k. 


iPrk^  +  _4/3) 


(ylPr)-4l3/ 


(36) 


where  Pr  =  i^/x  ,  the  Prandtl  number  (—0.7  for  air).  The  factor 
exp  [-iVt]  translates  the  resulting  inversion  so  that  is  replaced  by 
X3  -  Re  dZ3(0,  0)t  and  again  the  coordinate  system  offering  the  greatest 
utility  is  one  moving  upward  with  uniform  non-dimensional  velocity 
R.e  dZ3(0,  0)  =  Re  V3o(b/ff)3/2  .  jf  spectra  of  temperature  is  in¬ 
verted  and  expressed  in  this  system  the  result  is  (Appendix  C) 


8  =  8/n  V^Q  b 


3/2 


Re 


Pr 

yH 


l-4Pr/3  B 
y/Pr-4/3 


erf 


(y-1) 


/~y 

2,/~t+b 
^  Pr 


y/Pr-4/3  Bx, 


—  erf 
r 


^/4t/3  +b 


Performing  the  differentiation  yields 


(37) 


0 


Re 


Pr 
y  Hr 


erf 


(l-4Pr/3) 


where  T  =2,/ t+b  ,  T  =  2/4t73  +  b  and,  as  before, 
A  s  Pr  B 

S(u)  =  u  erf^  u  —  erf  u  . 
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If  we  let  b=  1/4,  the  length  scale  is  established  such  that  the  initial 
velocity  drops  by  a  factor  e  at  a  distance  from  the  center  of  the 
fireball. 


To  establish  the  velocity  scale  let 
\/  ,  X  small.  This  yields  0  =  2 


t  =  0,  r-0 
VjoRe  Pr 

yHT. 


and  note  that  erf  x-2x/ 
from  Equation  38  at 


r  -  0'^  for  b-1/4.  Define  6^^  as  the  value  of  6  at  t  =  0,  r  =  0  (at  the 
center  of  the  tireball)so  that 

2  v^Q  Re  Pr 

e  = - 

oc  yH 


Substituting  into  Equation  38  to  remove  the  Reynolds  number,  velocity 
product 


s=^e 

2  oc 


1  ,  r  3 

-  erf  - - T 

^  (y-4Pr/3)r 


{l-4Pr/3)S^^^ 


B, 


Along  the  h'>rizontal  midplane  x^  =  0  the  temperature  function  is 

P 


e  =  0  erf 

2p  oc 


2  v  2t  +  1/4 


where  we  have  taken  y  =  1 . 4,  Pr  =0.7  and  p  =  (x ^  +  x^  ) 


z^n 


This  function  is  relatively  constant  for  small  values  of  p/2v2t  +  l/4 
as  p  increases  indicating  the  presence  of  a  thin  "sheet"  of  high  tem¬ 
perature  air  in  the  horizontal  midplane  of  the  moving  coordinate 
system.  For  small  values  of  p/2  V  2t  +  1/4  the  temperature  for  x^  =  0  is 


0 


2=  v'2t  +  l/4  , 


giving  (using  e®=  T/T  ) 


\V8t+l 
oc 


^o  / 

Note  that  T^  is  the  temperature  of  the  atmosphere  far  from  the  fireball. 


^•'We  approach  the  moving  origin  along  a  horizontal  path. 
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SECTION  7 
DISCUSSION 


The  velocity  and  temperature  distributions  obtained  in  the  previous 
sections  result  from  an  extension  of  the  theory  of  turbulence  for  the 
final  period  of  decay.  The  extension  is  brought  into  the  theory  by  in¬ 
cluding  that  part  of  the  non-linear  terms  which  corresponds  to  coupling 
of  the  largest  eddies  with  eddies  of  all  other  sizes.  The  resulting  lin¬ 
earization,  which  is  carried  out  in  wave-vector  space,  has  the  effect 
of  translating  the  space  coordinates  in  the  prescription  of  the  velocity 
and  iemperature  in  a  w'ave-like  manner  depending  upon  the  velocity  of 
the  largest  eddies  in  the  initial  motion.  The  resulting  translational 
velocity  is  to  be  identified  with  the  rise  velocity  of  the  fireball  if  the 
initial  motion  has  only  an  upward  component.  A  convenient  coordinate 
system  for  describing  the  motion  and  temperature  in  the  case  of  more 
general  initial  conditions  is  one  which  is  moving  in  the  direction  of  the 
velocity  of  the  largest  eddies  in  the  initial  motion. 

From  a  physical  point  of  view  the  energy  source  most  effective  in  the 
motion  of  the  late-time  fireball  is  the  gravitational  field.  Its  origin 
is  in  the  downward  motion  of  those  parts  of  the  atmosphere  outside  the 
fireball.  This  may  be  shown  analytically  by  retaining  the  constant 
buoyant  force  term  through  the  transformation  to  wave-vector  space. 

The  term  in  the  transformed  equations  which  represents  the  buoyant 
force  pioduces  a  downward  velocity  which  is  effective  on  only  the 
largest  eddies  in  the  flow'  since  it  contains  the  Dirac  delta  function 
as  a  factor.  Attempts  to  continue  the  analysis  in  this  direction,  which 
is  taniamount  to  relaxing  the  assumption  of  pressure  equilibrium, 
were  unsuccessful  due  to  the  difficvilty  associated  with  the  interpreta¬ 
tion  of  integrals  w’hich  arose  containing  the  delta  function.  It  was  felt 
that  formal  application  of  the  rubies  fo-^  handling  the  delta  function  could 
lead  to  ambiguous  results  from  a  physical  po’nt  of  view'.  For  this  rea¬ 
son,  the  assumption  of  pressure  equilibrium,  which  states  that  the 
pressure  it  corresponding  altitudes  within  the  flow  is  equal  to  that  in 
an  undisturbed  exponential  atmosphere,  w'as  used  before  the  transforma¬ 
tion  was  applied  to  the  momentum  equation.  This,  in  effect,  removed 
the  energy  source  allowing  the  motion  to  proceed  in  accordance  with  the 
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theory  of  the  final  period  of  decay  of  turbulent  flow.  The  resulting 
effects  are: 

1.  To  remove  from  the  solutions  for  the  irrotational  part  of 
the  velocity  those  terms  which  correspond  to  propagating 
sound  waves  such  as  those  obtained  in  Appendix  A 

E.  To  remove  a  term  in  the  energy  equation  which  depends  on 
the  time  derivative  of  pressure 

3.  To  simplify  those  non-linear  terms  in  the  energy  equation 
which  involve  space  derivatives  of  the  pressure. 

Item  2  has  at  least  removed  the  effect  of  radiating  sound  waves  on 
the  temperature  distribution.  Item  3  has  an  effect  which  is  difficult 
to  determine  since,  in  e  S  S  x*C  6  >  it  is  a  linearization  process  brought 
about  by  the  fact  that  in  an  exponential  atmosphere  the  only  space 
derivative  of  the  logarithm  of  the  pressure  tTp  which  is  non-zero  is 
the  derivative  with  respect  to  the  vertical  direction  and  it  is  a  constant. 

The  overall  effect  on  the  motion  of  the  assumptions  considered  above 
is  to  smooth  out  irregularities  leading  to  a  deterministic  but  certainly 
approximate  description  of  the  velocity  and  temperature.  The  smooth¬ 
ing  which  occurs  prevents  this  method  from  providing  accurate  infor¬ 
mation  on  the  problem  of  radar  clutter  due  to  a  fireball.  The  method 
employed  does  provide  at  least  qualitative  information  on  the  velocity 
and  temperature  distribution.  An  extension  of  the  analysis  to  include 
the  influence  of  irregular  pressure  waves  maybe  required.  A  possible 
starting  point  is  to  employ  the  solutions  obtained  in  Appendix  A. 

The  form  of  the  solutions  obtained  in  this  report  have  suggested  an 
approach  to  the  problem  of  fireball  mixing  involving  the  use  of  Hermite 
functions.  Employing  a  Gaussian  distribution  of  upward  velocity  as  the 
initial  condition  the  solutions  found  are  expressed  in  terms  of  the  error 
function  and  its  derivatives.  These  derivatives  as  well  as  the  initial 
conditions,  may  be  written  as  Hermite  functions  which  have  several 
properties  which  may  be  useful  in  obtaining  solutions  of  the  fundamental 
equations.  For  example,  the  Fourier  transform  of  a  Hermite  function 
is  also  a  Hermite  function,  convolution  integrals  involving  Hermite  func¬ 
tions  may  be  performed  analytically,  a  series  of  Hermite  functions  is 

•  16 

capable  of  expressing  a  wide  range  of  generalized  functions. 

Another  alternative  approach  and  one  which  is  usually  employed  in  the 
general  theory  is  lo  determine  relationships  among  the  various  statistical 
features  of  the  dynamics,  such  as  correlation  coefficients,  but  to  include 
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better  approximations  to  the  governing  equations.  Some  of  these 
alternatives  are  considered  in  Appendix  B.  Since  no  direct  data 
exist  on  cor  ‘elation  coefficients  within  the  fireball,  the  spectra  of 
velocity  and  temperature  were  felt  to  be  more  useful  for  this 
application. 

Although  the  velocity  field  found  for  the  initial  conditions  used  here 
appears  to  represent  a  developing  vortex  ring,  the  temperature  in 
the  horizontal  midplane  of  the  moving  coordinate  system  has  a  form 
which  represents  a  vertically  thin  layer  or  membrane  of  high  temper¬ 
ature  air  within  the  ring  of  the  vortex  which  is  dropping  intemper- 
oture  due  to  ordinary  conduction.  This  occurs  because  the  velocity 
in  the  center  of  the  ring  is  small  in  the  moving  frame  of  reference 
and  the  thin  layer  rises  while  diffusion  of  energy  occurs  by  heat 
conduction.  At  other  positions  the  terms  involving  the  velocity  which 
appear  2.s  "source"  terms  in  the  energy  equation  act  to  increase,  by 
convective  transport,  the  temperature  outside  the  fireball  and  reduce 
it  on  the  inside.  In  this  way  the  fireball  "grows"  by  outward  move¬ 
ment  of  heated  air  which  mixes  with  the  atmosphere.  Only  part  of 
this  convective  effect  is  included  in  the  solutions  for  the  temperature 
found  in  this  report.  The  non-linear  term  Vj  SO/Sx^  which  -represents 
this  convective  effect  is  only  partially  included.  The  other  source 
term  represents  the  effect  of  expansion  and  compression  on  the  tem¬ 
perature  and  probably  is  less  important. 

If  the  solution  for  the  velocity  is  simply  substituted  into  the  energy 
equation  or  the  continuity  equation  (see  Footnote  page  13)  integrals 
are  found  which  represent  singular  solutions.  This  is  probably  re¬ 
lated  physically  to  the  impossibility  of  maintaining  exact  pressure 
equilibrium  in  the  flow.  The  full  understanding  of  this  and  many 
other  features  of  the  problem  must  be  reserved  for  later  investigation. 
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SECTION  8 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  conclusions  obtained  as  the  result  of  this  studv  are; 

1.  The  most  directly  useful  results  for  the  fireball  mixing 
problem  will  be  obtained  from  the  velocity,  temperature,  pres¬ 
sure,  and  density  rather  than  from  solutions  describing  cor¬ 
relation  cc-  ificients  involving  these  variables. 

2.  An  extension  of  the  theory  of  turbulence  in  the  final  period 
of  decay  may  be  obtained  by  including  that  part  of  the  non-linear 
inertia  terms  which  provides  a  translational  effect.  The  largest 
eddies  in  the  initial  motion  cc-nvect,  without  distortion,  the 
eddies  of  smaller  .?ize. 

3.  The  velocity  di:tribution  of  a  heated  mass  of  air  rising  in 
the  atmosphere  obtained  from  the  equations  of  change  of  a 
viscous,  heat  conducting,  compressible  gas  show  features  which 
correspond  in  several  respects  to  the  actual  characteristics  of 
the  late-time  fireball  of  a  nuclear  detonation. 

Recommendations  for  further  work  are: 

1.  A  more  detai’jd  comparison  of  the  results  obtained  in  this 
report  with  experimental  data  should  be  carried  out.  Different 
kinds  of  initial  conditions  should  be  employed. 

2.  Use  of  machine  computation  to  obtain  solut'  is  the  equations 
of  change  in  spectral  form  should  be  carried  cut.  This  approach 
would  begin  with  an  analysis  of  the  one-dimencional  case. 

3.  A  study  of  the  general  theory  of  turbulent  flow  as  related  to 
the  problem  of  the  late-time  fireball  should  be  undertaken.  This 
study  should  include: 

a.  Attempts  to  employ  the  solutions  of  Appendix  B  as  a 
starting  point  to  obtain  information  useful  for  radar  clutter 
problems. 
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b.  A  study  of  methods  of  generalizing  the  linearization 
process  employed  in  this  report. 

c.  An  investigation  of  the  use  of  higher  transcendental 
functions,  such  as  Hermite  functions  using  the  theory  of 

■n  •  iS 

Pansions . 
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APPENDIX  A 

ASYMPTOTIC  SOLUTIONS  OF  THE  LINEAR  EQUATIONS 


In  nou- dimensional  form,  after  Fourier  transformation  with  respect 
to  the  independent  space  variables,  the  equations  of  Moyal  are  of 
the  form 


+  dzf^  =  Re  dc!^^  (Z) 


+  I  dr.j^  +  Re^  ^  (dS+dQ)  =  Re  dCj^  (Z) 


+  i  k  Re  dZ,  =  Re  dC(S) 
at  k 


+  i  k  (y-  1)  Re  dZ,  +  V  ~  Q  =  Re  dC  (Q) 
at  K  1/ 


As  such,  as  is  usually  the  case,  successive  approximation  procedures 
are  suggested  for  the  two  limiting  cases  of  vanishing  Reynolds  number 
and  for  large  Reynolds  number.  Moyal  himself  wrote  down  a  solution 
for  the  homogeneous  case  neglecting  the  non-linear  convolution  integrals 
which  he  claimed  was  accurate  to  terms  of  order  1/  ..  In  dimensionless 
form  his  solutions  for  a  zero  mean  velocity  assume  the  form 
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We  have  noticed,  however,  that  this  approximation  has  the  following 
remarkable  property;  namely,  it  is  an  approximation  which  gives  the 
correct  sum  and  triple  product  of  the  three  roots  of  the  secular  equa¬ 
tion  for  the  last  three  coupled  equations  involving  the  longitudinal 
components,  but  it  does  not  correctly  give  the  sum  of  the  pair-wise 
products  of  the  roots.  As  such,  however,  it  would  form  a  good  basis 
for  computing  a  more  accurate  approximation  to  each  root  by,  for 
exaiTiple,  Newton's  method.  Since  we  were,  however,  interested  in 
trying  to  understand  at  least  some  of  the  simpler  aspects  of  the  non¬ 
linear  phenomena,  we  did  not  pursue  this  matter  further. 
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The  non-dimensional  form  of  these  equations  does,  however,  also 
suggest  a  successive  approximation  procedure  w'hich  in  principle 
would  yield  information  for  large  Reynolds  numbers  by  standard 
techniques.  It  is  to  be  noticed  that  the  first  equation  does  not  involve 
the  Reynolds  number  on  the  left-hand  side  while  the  coupled  set  are 
in  a  form  to  which  methods  have  been  developed  for  finding  a  solution 
to  the  homogeneous  equations  as  they  are,  for  example,  described 
in  Chapter  6  of  The  Theory  of  Ordinary  Differential  Equations,  by 
Coddington  and  Levison,  McGraw-Hill,  1955.  Once  such  a  funda¬ 
mental  set  of  solutions  were  known  along  with  the  trivial  solutions 
for  the  homogeneous  uncoupled  transverse  equations,  the  convolution 
integrals  on  the  right-hand  side  of  Moyal'  s  equations  could  be  evaluated. 
The  right-hand  side  of  the  equations  could  then  be  thought  of  as  known 
and  a  first  approximation  to  the  non-linear  effects  found  as  follows; 

The  uncoupled  equations  for  the  transverse  components  could  be  solved 
directly  by  the  method  of  variation  of  parameters  and  added  co  the  trivial 
general  solution  to  the  homogeneous  part  of  these  equations.  If  the  right- 
hand  side  of  the  coupled  equations  describing  the  component  along  k  is 
denoted  by  the  vector  b(t,  k)  and  if  the  fundamental  matrix  for  a  large 
value  of  the  Reynolds  number,  found  by  solving  the  homog  neous  coupled 
equations  by  a  procedure  we  will  outline,  were  denoted  by  $(t,k)  then  the 
vector  function  0(t,k)  defined  by  the  relation  [Equation  3.  1  page  74  of 
Coddington  and  Levinson] 


^(t,k)=  $(t,k) 


(s,  k)  b(s ,  k)  ds 


wt  id  be  a  solution  of  the  non-homogeneous  equations  vanishing  for  t^  . 
To  this  could  be  added  the  general  solution  of  the  homogeneous  equations 
already  found  in  order  to  have  a  first-order  approximation  to  the  non¬ 
linear  effects. 


Unfortunately,  time  did  not  permit  us  to  explore  this  mode  of  attack  in 
detail  although  we  did  make  a  beginning  in  this  direction.  Since  there 
is  some  possibility  that  work  in  this  area  will  be  funded  on  a  long-range 
basis  in  the  future,  it  appears  to  us  to  be  worthwhile  to  record  the 
progress  obtained  to  date  by  this  approach.  If  for  simplicity  we  let 

=  Kfv  -  ,  where  Pr  is  the  Prandtl  number,  then  the  coupled 

homogeneous  equations  of  Moyal  can  be  written  in  the  matrix  form 
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That  is  in  the  general  matrix  form 


d  X 
dt 


Re 


>0  + 


1 

Re 


X 


If  Bq  is  a  non- singular  matrix  such  that  where 

is  a  diagonal  matrix  with  distinct  roots  that  can  be  found,  as  we  will 
show  to  be  the  case  for  Moyal's  equations,  then  the  transformation 
X  =  BqY  leads  to  a  system  of  the  form 


Re 


1 


(A-l) 


For  Moyal's  equations  the  situation  is  particularly  simple  for  it  is  pos¬ 
sible  to  find  a  matrix  Bq  with  the  above  properties  which  is,  in  fact, 

(')  -1 

independent  of  the  time.  In  this  situation  A  j  =  B  Aj  B^  .  For  sim¬ 
plicity  we  will  again  denote  this  by  Aj  . 


We  will  now  briefly  describe  the  method  used  to  obtain  a  formal  solu¬ 
tion  to  Equation  A-l.  This  description  is,  in  reality,  a  specialization 
of  that  given  in  Chapter  6  of  Coddington  and  Levinson,  but  it  seems 
worthwhile  to  include  it,  for  the  case  at  hand  is  considerably  simpler 
than  the  general  case  treated  by  them.  We  seek  a  solution  for  ^  in  the 
form 

Jr'  =  P  exp  Q  ,  (A -2) 

where  P(t)  and  Q(t)  are  matrices,  Q(t)  being  assumed  diagonal.  For 
P(t)  we  assume  an  expansion  of  the  form 

P(t)  =  Po(t)  V  P.(t)  +  P2(t)  +  .  .  .  . 

Re‘' 

while  for  Equation  A-l  it  suffices  to  take  Q  =  Re  Ooit)  +  QjCt)  where, 

as  we  shall  show,  O  (t)  =  D„,  a  constant  matrix  and  where  Qi  is  also 

o  o  ^ 

a  diagonal  matrix  by  assumption. 
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Performing  the  indicate  differentiation  and  substituting  into  Lcnation 
A-1  we  obtain  after  canceJlation  of  the  exponential  factor 


P^  + 

Pi  +—  Pp 
^  „  2  2 

+ 

1  p  +  _i_ 

p,  +  — ,  p,  +  ...) 

®  Re 

1  o  Re 

Re 

Re 

Re  Qq  ^  Gjl  -  Re  I  I  D„  + 


f  1  1 

■P.  +  P ,  +  - -  +  .  .  . 

Re 


Re  1)  1  o  Re  ‘  1 


Rearranging  according  to  decreasing  powers  of  the  Reynolds  number, 
the  above  becomes 


PoQ;  -  Do  Po  +  (R'l  [p;  +  Pi  PoQ]'  -  DoPl  '  ^iPol 

!Pl'"p^Q;"Pl  q; -DoPz-AiP,) 

Pz  ^3  ^o  ^2  ^2  -  ^0^3  “"^1^2^  -f  ,  .  .  =  0  . 


Re 


By  equating  the  various  powers  of  the  Reynolds  number  to  zero  we 
obtain  the  following  infinite  set  of  matrix  equations: 


P„Q  -DP 
o  o  o  o 


0 


p;  +  Pi°; "  Poq;  -  doPi-a,Po  =  o 
P;  "PjQ^^PiQ;  -  DoPj-AiP,  =  0 
Pz  +P3°o  +  PeQi'  -  DoPz-AjPj  =  0 


The  first  of  these  is  easily  satisfied  by  taking  P^  =  I  =  the  identity 
matrix  which,  of  course,  implies  that  P^  =0  and  setting  0^=0  . 

To  satisfy  the  second  equation  it  is  first  rewritten  in  the' form 

in  which  use  has  been  made  of  the  information  already  obtained.  Now 
it  is  clear  that  there  are  no  diagonal  terms  on  the  left-hand  side  of 
the  above  equation.  Thus,  the  diagonal  terms  on  the  right-hand  side 
must  cancel  so  that  if  we  let 

Ai  =  (aij) 
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we  See  that  since  Qj  is  diagonal,  by  assumption  Qj  =  {ajj  S^j)  where 
6ij  is  the  Kronecker  delta  equal  to  one  for  equal  indices  and  equal  to 
zero  for  different  indices.  For  i#  j  we  then  find  that  if  Pj  (p^j), 
i  ^  j  then 

Pij  =  ■ 


where  ,  Xj  are  distinct  roots  of  =  (X^  6ij)-  In  order  to  determine 
the  diagonal  terms  of  Pj  ,  Pj  is  replaced  by  Pj  =  Pj+Dj  where  Dj 
is  as  vet  an  unknown  diagonal  matrix  and  where  P,  is  assumed  to  have 
a  zero  diagonal.  P'j'  is  clearly  still  a  solution  of  the  second  equation 
along  with  Pj  for  such  an  arbitrary  Dj  .  To  determine  Dj  we  use  the 
next  equation  rev/ritten  in  the  form 


o ' 
•  1 


+  d'  +  PjOj  -  A 


1^1 


=  A,D, 


-D,Ql'  +P2D 


o 


-DP 

o 


7 


In  this  form  there  are  no  diagonal  terms  on  the  right-hand  side  and 
since  Pj  had  a  zero  diagonal,  it  follows  that  the  diagonal  matrix  Dj 
must  be  given  by  the  diagonal  terms  of  AjPj  -  PjQ]  •  This  same 
process  can  be  continued  to  higher  order  but  we  shall  not  carry  it  out. 
We  w'ill  continue  our  comments  on  the  general  theory  with  the  observa¬ 
tion  that  for  the  equations  of  Moyal  ,  A  j  hence  Oq  ,  Pj  ,  Qj  and 
Dj  all  turn  out  to  be  matrices  independent  of  the  time.  This  means 
that  the  above  procedure  will  produce  in  this  particular  case  a  partic¬ 
ular  simple  type  solution  of  the  form 

1  =  [I  +  {Pi  t)]  1  {^i  ^ij  ^  {siij  ^ij)  1^  ’ 


or 


x  =  exp  jRe  (Xj  6ij)  +  (a^j  6ij)t| 


+  —  (B^Pi  +  B^Di  t)  exp  [Re  (Xi  6^3)  +  {a-  6ij)jt 


To  apply  this  theorm  to  the  equation  of  Moyal  involving  the  component 
along  k,  it  is  first  necessary  to  determine  the  eigenvalues  of  Aq  from 
the  equation 

det  (XI  -  Aq)  =  0 
or 
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=  0 


A  simple  computation  yields  the  three  roots 


^o  =  0 


a 


—  i  a  k 


X^  =  i  a  k 


z  2  z  z 

We  note  the  two  relations  Xj  =  ^2  =  -  a  k  ,  iXj  -  1X2  =  2  a  . 

In  order  to  determine  the  matrix  with  the  property  that 

=  D 

o  O  o  o 

we  use  the  remark  on  page  176  of  Coddington  and  Levinson  and  take 
for  the  k-th  column  of  Bq  a  multiple  of  the  cofactors  of  a  row  of 
Ap  -  Xj^I  .  Specifically,  we  chose  for  Bp  the  matrix 


i  X 


1 


i  X- 


-Ik  k 

1  k(v-l)  k(y-lk 


We  next  verify  the  relation  A^  B^  =  Ep  Dp  by  direct  computation, 
making  use  of  the  simple  properties  possessed  by  the  characteristic 
roots  noted  previously.  Thus 
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These  relations  make  explicit  the  proper  form  for  Bq  and  for 

both  right  and  left  matrix  multiplication.  We  next  complete 

=  Bq^  AjBq  in  two  successive  steps. 
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We  next  note  that  Bq^  can  be  written  as 
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Therefore 
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Qi 


and 


Q 

e 


ReQo+Qj 
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where  A,  B,  C  are  constants 
of  integration. 


If  we  now  define  e^  as  the  logarithm  of  a  matrix  J  (see  Coddi  <?ton 
and  Levinson  page  65),  i.e.  , 


and  think  of  Q  in  the  form 
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log  Ml  0  0 

0  log  M2  0 

0  0  log  M  3 
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where  M  =  -  —  k - — 

^  Z 

Thus,  if  we  keep  only  the  first  term  in  the  solution 
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we  find  that 
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(A-3) 


(A-4) 


(A-5) 


dZ^  =  aktB-C) 
dS  I  =  -A+k  (B+C) 
dQ  1^  =  A  +k(y-l}  (B+  C) 

Add  Equations  A- 7  and  A -8 

*  I 

dP  io  =  yk  (B  T  C  ) 
and  with 

'^'^klo=  ak(B  -  C)  . 


(A-6) 
{A -7) 
(A-8) 
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This  yields  for  A,  B,  and  C 

dPlo/yk  +  dZj^l  jaik-  ZB 

dP|o/yk-dZ^lJak  =  2C 

A  =  dQ  i  -  (v- 1)  dP  iy  . 
lo  o/ 


Let  a  =  r  I  +-|  (y-l) 


Substituting  back  into  the  set  of  Equations  A-1  through  A- 3, 

(-ia  kRe  -ak^  )t 


dZk=ciZj^„e 


dS  = 


-dQ  +  dP 

o  y  o 


(—  i  a  k  Re  -  ak^)t 


-0k^t 
e  + 


/dZ 


k  dP 
o  .  o 


2a 


2v 


dP 


dZ, 


27 


2a 


1  (+i  a  k  Re  -  Qfk^)t 

e 


dQ  =  (dQ„-^  dP  1  +(y-l) 


dZ.  dP 
^  lo  ^  o 


(-iakRe  -  ak'^)t 


(■KakRe-ak^)t 


+  (r-i) 


dP 


2y 


2a 

dZj^l  , 


2y 


2a 
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APPENDIX  B 

REVIEW  OF  THE  LITERATURE 


The  mathematical  theory  of  turbulence  has  a  long  history.  Except  in 
the  simplest  circumstances  detailed  understanding  is  far  from  complete. 
So  far,  most  success  which  correlates  with  experimental  results  has 
been  achieved  for  the  special  case  of  homogeneous  isotropic  incompress¬ 
ible  turbulence  as  it  is,  for  example,  generated  far  down  stream  from 
wire  grid  in  a  wind  tunnel.  The  theory  for  this  case  is  based  on  the 
appropriate  Navier-Stokes  equations  from  which  various  correlation 
functions  are  constructed  by  taking  averages  on  the  equations  them¬ 
selves  rather  than  on  their  solutions  which  are  unknown.  This  process, 
unfortunately,  always  leads  to  an  open  set  of  moment  equations  because 
of  the  presence  of  the  non-linear  terms.  Thus,  with  this  approach 
some  additional  hypothesis  is  necessary  in  order  to  close  the  system. 

In  the  simplest  case  this  merely  consists  in  neglecting  all  correlation 
functions  of  higher  than  the  second  order;  but  a  great  many  more  subtle 
approaches  have  also  been  used,  including  the  quite  widely  used  quasi- 
Gaussian  approximation  which  relates  fourth  order  correlation  functions 
to  those  of  second  order.  The  use  of  spectral  analysis  via  methods  of 
Fourier  transformation  have  also  been  extensively  exploited  and  much 
work  has  been  devoted  to  the  hypotheses  introduced  by  Heisenberg  and 
Kolmogoroff  (the  so-called  locally  homogeneous  and  locally  isotropic 
turbulence),  whose  purpose  was  to  determine  the  decay  laws.  We  will 
not  enter  into  any  of  these  matters  here  in  detail  for  they  are  nicely  de¬ 
scribed  in  the  book  by  G.  K,  Batchelor  (The  Theory  of  Homogeneous 
Turbulence,  Cambridge  Univ,  Press,  1953)  as  well  as  in  the  book  by 
L.  Agostini  and  J,  Bass  (Les  Theories  de  !a  Turbulence,  Publ.  Sci. 
Tech,  Ministere  L'air,  no.  237,  1950). 

In  addition  to  the  approach  described  above  there  are  at  least  two  other 
distinct  methods  of  attack  which  are  conceptually  more  appealing.  The 
first  consists  in  the  deduction  of  a  functional  equation,  equivalent  to 
an  infinite  set  of  equations,  whose  solution  where  it  is  k.nown  would 
yield  exact  expressions  for  the  various  correlation  functions  of  interest. 
This  approach  was  initiated  by  E,  Hopf  (Statistical  Hydromechanics 
and  Functional  Calculus,  J.  Rati.  Mech.  Anal.  ,  vol  1,  1952,  pp  87-123) 
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in  a  space  formalism  aad  recently  extended  to  a  space  time  formalism 
by  R.  Lewis  and  R.  Kraichnan  (A  Space-Time  Fvmctional  Formalism 
for  Turbulence,  Comm,  Pure  and  App.  Math.  15,  1963,  p  363). 

The  second  method  involves  the  representation  of  the  various  quantities 
such  as  velocities,  density,  temperature,  etc,  ,  in  terms  of  Fourier 
transforms  of  singular  random  functions  or  fields.  The  basic  concepts 
are  briefly  discussed  in  Batchelor's  book  on  pages  28-33  and  in  the 
basic  paper  by  Moyal  on  which  most  of  this  report  is  based  (The 
Spectra  of  Turbulence  in  a  Compressible  Fluid;  Eddy  Turbulence, 
and  Random  Noise,  Proc,  Camb.  Phil,  Soc.  48,  1952,  page  329).  The 
relevant  mathematical  ba.ckground  is  discussed  by  Moyal  (Stochastic 
Processes  and  Statistical  Physics,  J,  Royal  Statistical  Soc.  ,  B,  11, 
1949,  p  150).  This  approach  has  also  been  the  subject  of  much  Russian 
work  as  may  be  seen  by  reference  to  the  article  by  A.  Yaglom  (Some 
Classes  of  Random  Field  in  n-Dimensional  Space,  Related  to  Stationary 
Random  Processes,  Theory  of  Probability  and  its  Applications,  2,  no. 

3,  1957,  p  273),  While  this  work  too  has  largely  been  confined  to 
characterizations  involving  incompressible  flows,  the  paper  by  Moyal 
is  one  of  the  few  that  we  have  discovered  which  lays  a  foundation  on 
which  one  could  base  a  long-range  systematic  attack  on  the  problems 
of  interest.  In  fact,  our  search  of  the  literature^  for  material  relat¬ 
ing  to  compressible  turbulent  effects  proved  most  disappointing,  al¬ 
though  not  unexpectedly  so,  because  of  the  intrinsic  difficulty  of  the 
problem.  In  particular,  Moyal  achieves  a  partial  separation  in  wave 
vector  space  of  the  eddy  forming  part  of  the  velocity  (the  transverse 
part  of  the  singular  random  field)  from  the  remainder  (the  longitudinal 
part)  which  allows  a  somewhat  easier  physical  interpretation  of  the 
involved  complicated  phenomena.  Before  entering  into  Moyal's  re¬ 
sults,  it  is  perhaps  worthwhile  to  review  in  brief  some  of  the  facts 
concerning  the  nature  of  turbulence  as  expressed  by  three  distinguished 
mathematicians:  J,  Leray,  E.  Hopf,  and  J.  von  Neumann. 

We  begin  this  review  by  quoting  in  free  translation  some  excerpts 
from  Leray's  paper  (Sur  le  Mouvement  d’un  Liquide  Visqueux  Emplis- 
sant  L'espace,  Acta.  Math.  ,  63,  1934,  p  193). 


*A  rather  thorough  and  systematic  literature  search  was  performed 
using  Mathematical  Reviews,  Science  Abstracts,  and  DDC. 
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"The  theory  of  viscosity  of  liquids  is  governed  by  the  equa¬ 
tions  of  Navier=  It  is  necessary  to  justify  a  posteriori 
this  hypothesis  by  establishing  the  existence  theorem  which 
follows.  There  exists  a  solution  of  the  equations  of  Navier 
which  corresponds  to  a  given  arbitrary  initial  set  of  values 
of  the  velocities.  Oseen  tried  to  prove  this  —  he  only  suc¬ 
ceeded  in  establishing  the  existence  of  such  a  solution  for 
a  period  of  time,  perhaps  very  short,  succeeding  the  initial 
instant.  Now  one  can  verify  that  the  total  kinetic  energy 
of  the  fluid  remains  bounded,  but  it  does  not  seem  possible 
to  deduce  from  this  the  fact  that  the  flow  remains  regular 

"It  is  not  paradoxial  to  suppose  that  the  cause  regulating 
the  flow  —  the  dissipation  of  energy — is  not  sufficient  to 
keep  the  second  derivatives  of  the  velocity  components 
bounded  and  continuous  with  respect  to  their  coordinates; 
now  the  theory  of  Navier  assumes  that  the  second  deriva¬ 
tives  are  bounded  and  continuous,  Oseen  himself  comment¬ 
ed  on  the  unnatural  character  of  this  hypothesis;  at  the 
same  time  he  showed  how  to  express  the  laws  of  fluid  me¬ 
chanics  in  terms  of  integrodifferential  equations  in  which 
only  the  components  of  velocity  and  their  first  derivatives 
with  respect  to  the  coordinates  occur.  In  tue  course  of 
this  work  I  will  consider  a  tern  of  relations  which  is 
equivalent  to  the  integrodifferential  equation  of  Oseen,  com¬ 
pleted  by  an  inequality  expressing  the  dissipation  of  energy. 
These  relations  are  deduced  from  the  equations  of  Navier 
by  the  aid  of  integration  by  parts  which  makes  the  highest 
order  derivatives  disappear.  (Today  these  are  the  so- 
called  weak  solutions.  And  if  I  have  not  succeeded  in 
establishing  the  existence  theorem  stated  above,  I  have 
nevertheless  proved  the  following:  the  relations  in  ques¬ 
tion  (i.  e.  ,  the  weak  form)  always  possess  at  least  one 
solution  which  is  defined  for  an  unlimited  time  and  which 
corresponds  to  a  given  initial  state  of  velocity.  Perhaps 
this  solution  is  not  regular  enough  to  possess  second  deriv¬ 
atives  which  are  bounded  at  each  instant  of  time;  nonethe¬ 
less  it  is,  in  a  proper  sense  of  the  term,  a  solution  of  the 
equations  of  Navier;  I  propose  to  call  them  turbulent  solutions. 

"If  I  had  succeeded  in  constructing  the  solutions  of  the 
Navier  equations  which  became  irregular,  1  would  then 


*  Authors*  note. 
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be  able  to  affirm  that  turbulent  solutions  exist  which  are 
not  at  all  simply  reducible  to  regular  solutions.  But  if 
this  is  false,  the  concept  of  a  turbulent  solution  will  not 
lose  its  interest  since  it  is  bound  to  play  a  large  role  in 
the  study  of  viscous  liquids;  it  is  well  designed  to  handle 
the  problem  of  physics  and  mathematics  where  the  physical 
causes  of  regularity  are  not  sufficient  to  include  the  hypo¬ 
theses  of  regularity  which  are  consequences  of  the  deduced 
equation;  one  can  apply  the  considei ations  presented  here 
to  these  problems. 

"We  note  the  following  two  facts: 

"Nothing  permits  one  to  assert  the  uniqueness  of  the 
turbulent  solution  which  corresponds  to  a  given  initial 
state; 

'The  solution  which  corresponds  to  a  state  sufficiently  near 
to  a  steady  state  can  never  become  irregular.  " 

To  express  the  philosophy  of  E.  Hopf  we  quote  from  A  Mathematical 

Example  Displaying  Features  of  Turbulence,  Comm,  Pure  and  App. 

Math.  ,  1,  1948,  p  303. 

"It  is  convenient  to  visualize  the  solutions  in  the  phase  space 
of  the  problem.  A  phase  or  state  of  the  fluid  is  a  vector 
field  ...  in  the  fluid  space  that  satisfies  the  (Navier -Stokes) 
equations  and  the  boundary  conditions.  The  totality  Qof 
these  phases  is  therefore  a  functional  space  with  infinitely 
many  dimensions.  A  flow  of  the  fluid  represents  a  point 
motion  in  and  the  totality  of  these  phase  motions  form 
a  staticnarv  flo^v  in  the  phase  space  f?  ,  which,  of  course, 
is  to  be  distinguished  from  the  fluid  flow  itself.  What  is 
the  asymptotic  future  behavior  of  the  solutions,  how  does 
the  flow  behave  for  t  -•*  ,  and  how  does  this  behavior  cliange 
as  u  (the  viscosity)  decreases  more  and  more?  How  do 
the  solutions  which  represent  the  observed  turbulent  mo¬ 
tions  fit  into  the  phase  picture?  .  .  , 

"The  observational  facts  about  hydrodynamic  flow  reduced 
to  the  case  of  fixed  side  conditions  and  with  u  as  the  only 
variable  parameter  are  essentially  these;  for  U  sufficiently 
large,  1/  >  Uq  »  the  only  flow  observed  in  the  long  run  is  a 
stationary  one  (Laminar  flow).  This  flow  is  stable  against 
arbitrary  initial  disturbances.  Theoretically  the  correspond¬ 
ing  exact  solution  is  kno'"n  to  exist  for  every  value  of  V  >  0 
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and  its  stability  in  the  large  can  be  rigorously  proved, 
though  only  for  sufficient  large  values  of  U  .  The  cor¬ 
responding  flow  in  phase  space  Q  thus  possesses  an  ex¬ 
tremely  simple  structure.  The  laminar  solution  rep¬ 
resents  a  single  point  in  0  invariant  under  the  phase  flow. 

For  U>  every  phase  motion  tends  as  t-®  toward  this 
laminar  point.  For  sufficiently  small  values  of  U  ,  how¬ 
ever,  the  laminar  solution  is  never  observed.  The  tur¬ 
bulent  flow  observed  instead  displays  a  complicated  pat¬ 
tern  of  apparently  irregular  moving  'eddies'  of  varying 
sizes.  The  view  widely  held  at  present  is  that  for  u  >  0  , 
having  a  fixed  value,  there  is  a  'smallest  size'  of  eddies 
present  in  the  fluid  tending  to  zero  as  v-*0.  Thus,  mac- 
roscopically,  the  flow  has  the  appearance  of  an  intricate 
chance  movement  whereas,  if  observed  with  sufficient 
magnifying  power,  the  regularity  of  the  flow  would  never 
be  doubted. 

"The  qualitative  mathematical  picture  which  the  author 
conjectures  to  correspond  to  the  known  facts  about  hydro- 
dynamic  flow  is  this;  to  the  flows  observed  in  the  long  run 
after  the  influence  of  the  initial  conditions  has  died  down 
there  corresponds  certain  solutions  of  the  Navier-Stokes 
equations.  These  solutions  constitute  a  certain  manifold 
M  =  M(p)  in  phase  space  invariant  under  the  phase  flow. 
Presumably  owing  to  viscosity  M  has  a  finite  number  N  = 

N(p)  of  dimensions.  This  effect  of  viscosity  is  most  evident 
in  the  simplest  case  of  U  sufficiently  large.  In  this  M 
is  simply  a  single  point  N  =  0  .  Also  the  complete  stability 
of  M  is,  in  this  simplest  case,  obviously  due  to  viscosity. 

On  the  other  hand,  for  smaller  and  smaller  values  of  U  , 
the  increasing  change  characteristic  of  the  observed  flow 
suggests  that  N(P)  -*®  monotonically  as  iz-O.  This  can  hap¬ 
pen  only  if  at  certain  'critical'  values 

•  -0 

the  number  N(P)  jumps.  The  manifold  M(l/)  itself  presvimably 
changes  analytically  as  long  as  no  critical  value  is  passed. 
Now  we  believe  tha+  when  U  decreases  through  such  a  value 
,  a  continuous  branching  phenomenon  occurs.  The  mani¬ 
fold  M(P)  of  motions  observed  in  the  long  rvm  (more  precisely 
its  einalytical  continuation  for  V  <  loses  its  stability. 
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The  notion  of  stability  here  refers  to  the  whole  manifold 
and  not  to  the  single  motions  contained  in  it.  The  loss  of 
stability  implies  that  the  motions  on  the  analytically  continued 
M  are  no  longer  observed.  What  we  observe  after  passing 
is  not  the  analytical  continuation  of  the  previous  M  but 
a  new  manifold  M{v)  continuously  branching  away  from 

and  slightly  swelling  in  a  new  dimension.  This  new 
M(t'K)  takes  over  stability  from  the  old  one.  Stability  here 
means  that  the  'majority'  of  phase  motions  tend  for  t-*“  to¬ 
ward  M(t')  ,  We  must  expect  there  is  a  'minority'  of  ex¬ 
ceptional  motions  that  do  not  converge  toward  M  (for  in¬ 
stance  the  motions  on  the  analytical  continuation  of  the  old 
M  and  of  all  the  other  manifolds  left  from  all  the  previous 
branchings).  The  simplest  case  of  such  a  bifurcation  with 
corresponding  change  of  stability  is  the  branching  of  a  per¬ 
iodic  motion  from  a  stationary  one.  This  case  is  clearly 
observed  in  the  flow  around  an  obstacle  (transition  from 
the  laminar  flow  to  a  periodic  one  with  periodic  discharge 
of  eddies  from  the  boundary).  The  next  simplest  case  is 
the  branching  of  a  one-parameter  family  of  almost  periodic 
solutions  from  a  periodic  one,  " 

In  this  same  paper,  Hopf  gave  a  simple  example  of  a  mathematical 
model  which  displayed  the  features  described  above.  Later  in  a  paper 
entitled,  On  Non-Linear  Partial  Differential  Equations,  which  appeared 
in  the  Lecture  Series  of  the  Symposium  on  Partial  Differential  Equa¬ 
tions,  published  by  the  University  of  Kansas  in  1957,  he  gave  a  still 
simpler  example.  However,  suggestive  as  these  examples  are,  in 
no  sense  are  they  complicated  enough  to  lead  to  any  real  insight  into 
the  nature  of  actual  turbulent  flow  nor  do  they  suggest  any  attack  on 
the  main  problem  even  for  the  incompressible  case. 

Finally,  we  will  quote  a  few  excerpts  from  a  privately  circulated  re¬ 
port  by  the  late  John  von  Neumann,  dated  IZ  June  1950.  This  report 
was  made  available  to  us  through  the  courtesy  of  Professor  G,  Uhlenbeck 
of  the  Rockefeller  Institute, 

"It  is,  by  the  way,  obvious  that  one  has  to  talk  not  of  one 
turbulent  solution,  but  of  many  turbulent  solutions.  The 
problems  which  lead  to  turbulence  are  usually  such  that 
the  outer  conditions  of  the  problem  do  not  change  with  time, 
whereas  the  turbulent  solution  itself  invariably  does.  A 
translation  of  the  whole  solution  along  the  time- axis  is 
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therefore  always  possible  and  produces  infinitely  many 
solutions  of  the  same  problem,  all  possessing  the  same 
turbulent  characteristics  as  the  original  solution.  Indeed, 
there  is  every  reason  to  believe  that  in  addition  to  this, 
many  other  principles  of  transforming  these  solutions  exist. 
There  is  probably  no  such  thing  as  a  most  favored  or  most 
relevant,  turbulent  solution.  Instead,  the  turbulent  solu¬ 
tions  represent  an  ensemble  of  statistical  properties,  which 
they  all  share,  and  which  alone  constitute  the  essential  and 
physically  reproducible  traits  of  turbulence. 

"In  the  classical  setup,  statistical  mechanisms,  the  system 
under  consideration  is  isolated  from  the  outside  world. 

Its  energy  is  fixed,  and  one  is  concerned  with  the  to-and-fro 
energy  exchanges  between  degress  of  freedom.  The  clas¬ 
sical  theory  allows  to  infer  from  this  that  equipartition  of 
energy  between  all  degrees  of  freedom  prevails.  In  the 
present  setup,  on  the  other  hand,  the  system  is  'open'  at 
both  ends,  energy  is  being  supplied  as  well  as  dissipated. 
The  two  'ends'  do  not,  however,  lie  in  ordinary  space,  but 
in  its  Fourier  transform.  More  specifically:  The  supply 
of  energy  occurs  at  the  macroscopic  end — it  originates  in 
the  forced  motions  of  macroscopic  (bounding)  bodies,  or 
in  the  forced  maintenance  of  (again  macroscopic)  pressure 
gradients.  The  dissipation,  on  the  other  hand,  occurs 
mainly  at  the  microscopic  end,  since  it  is  ultimately  due 
to  molecular  friction,  and  this  is  most  effective  in  flow- 
patterns  with  high  velocity  gradients,  that  is,  in  small 
eddies.  (While  the  main  dissipative  flow- patterns  are 
microscopic,  they  are  not  necessarily  molecular,  .  .  .  .  ) 
That  is:  A  flow  of  energy  is  taking  place  from  the  sources 
which  are  situated  at  the  low  frequencies  to  the  sinks  that 
lie  in  the  high  frequencies.  (The  frequencies  are  to  be 
understood  in  space,  and  not  in  time,  )  Thus,  the  statistical 
aspect  of  turbulence  is  essentially  that  of  a  transport  phe¬ 
nomenon  (of  energy) — transport  in  the  Fourier-transform 
space. 

"Having  understood  this  decisive  trait — that  turbulence  is 
not  a  matter  of  ergodic  distribution  of  a  fixed  amount  of 
energy,  but  the  transport  of  a  fixed  flow  of  energy  from 
sources  in  the  low  frequencies  to  sinks  in  the  high  fre¬ 
quencies  in  the  Fourier-transform  space — it  is  now  pos¬ 
sible  to  take  the  next  step. 
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"This  step  consists  of  noting  th.  t  generally  valid  and  sim¬ 
ple  phenomena  should  be  expected  in  the  domain  of  pure 
flow  (energy  flow,  in  the  Fourier -transform  space),  re¬ 
mote  from  both  the  sources  and  the  sinks,  that  is,  in  the 
intermediate  frequencies.  Indeed,  the  sources,  as  well 
as  the  sinks  cause  specific  complications,  depending  on 
more  complicating  factors  with  which  it  is  advisable  to 
deal  separately  and  subsequently.  A  simple  and  unique 
phenomenon,  which  is  of  the  same  nature  under  all  pos¬ 
sible  conditions  of  turbulence,  is  to  be  expected  in  the 
domain  of  pure  energy  flow  only,  as  indicated  above. 

"This  circumstance  was  recognized  and  primarily  dealt 
with  in  the  work  of  group  C.*  The  discussion  which  fol- 
lov/s  will  accordingly  begin  with  emphasizing  this  aspect. 

The  variant  to  be  presented  is  a  composite  one,  using 
elements  of  the  approaches  of  several  authors. 

"The  intermediate  (spatial)  frequency  range  corresponds 
clearly  to  an  intermediate  size  range,  that  is,  to  an  inter¬ 
mediate  vorticity- and -eddy -size  range.  Note  again:  This 
range  has  to  be  defined  as  being  too  small  to  be  affected  by 
the  bounding  solids  (energy  sources),  but  at  the  same  time 
too  large  to  be  affected  by  high-  (molecularly  viscous)  dis¬ 
sipation  eddies  (energy  sinks).  Hence  there  should  be  no 
influence  of  linear  dimensions  (of  the  macroscopic  system): 

L  ,  and  also  no  influence  of  the  (molecular)  viscosity;  u  . 

Note  further:  Yet,  these  enter  into  the  Reynolds  number 

R  =  ,  which  was  heretofore  viewed  as  the  decisive  quan¬ 

tity:  Thus  the  only  quantity  that  can  matter  is  the  rate  of 
dissipation  W  (  =  energy/unit  mass  X  unit  time). 

"The  following  circumstance  is  also  worth  noting:  It  need 
not  be  true  that  the  mere  cessation  of  stability  of  the  lam¬ 
inar  flow  is  already  equivalent  to  the  onset  of  a  properly 
turbulent  regime.  It  is  quite  normal  to  expect  that  for 
Reynolds  numbers  only  slightly  in  excess  of  the  stability 
limit  only  one  unstable  mode  of  perturbation  exists.  The 
non-laminar  form  of  motion  which  then  develops,  will  there¬ 
fore  still  be  one  of  relative  simplicity  and  not  at  all  resem¬ 
bling  the  familiar,  highly  involved,  pattern  of  turbulence. 

HeAutiior'::  note:  The  work  of  A,  N.  Kolmogoroff,  L.  Onsager,  C.  F. 

von  V/eizrsricker  a  id  W.  Heisenberg. 
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(All  available  conclusive  results  in  laminar- perturbation 
stability  theory  confirm  this  .  .  .  .  )  The  development  of 
such  an  unstable  mode  of  perturbation,  that  is,  of  one  which 
increases  with  time,  nevertheless  may  start  a  chain  of 
events  which  ultimately  leads  to  turbulence.  The  plausible 
picture  of  this  evolution  may  be  somewhat  like  this:  As 
long  as  the  (unstable)  perturbing  mode  has  not  developed 
to  a  large  size,  the  stability  properties  of  the  perturbed 
(that  is,  distorted)  laminar  flow  will  be  essentially  the  same 
ones  as  those  of  the  unperturbed  (original)  one.  That  is, 
no  new  unstable  modes  come  into  play.  If  the  perturbing 
mode  reaches,  however,  a  certain  finite  size  —  at  which 
time,  of  course,  the  applicability  of  linear  theory  will  have 
ceased  and  a  non-linear  discussion  will  have  become  neces¬ 
sary —  the  stability  properties  may  shift.  The  now  signif¬ 
icantly  perturbed  (distorted)  flow  may  now  acquire  further 
unstable  modes  which  will  also  begin  to  grow.  They  are 
likely  to  catch  up  with  the  primary  unstable  mode  in  size 
and  importance,  because  the  primary  mode  is  already  in 
the  non-linear  range  and  therefore  probably  no  longer  in¬ 
creasing,  while  the  secondary  modes  will  be  exponentially 
increasing  as  long  as  they  are  small,  that  is,  until  they 
themselves  acquire  a  finite  size  and  leave  the  linear  ramge. 
As  long  as  '■hese  secondary  modes  are  smalU  they  will  cause 
no  change  of  the  stability  properties.  When  they  become 
large,  however,  they  may  in  their  turn  modify  the  stability 
situation  and  give  rise  to  ternary  unstable  modes,  etc.  It 
is  clear  how  this  mechanism  successively  produces  the 
degree  of  involvement  which  characterizes  turbulence  in 
the  ordinai^y  sense.  It  is  plausible  to  say  that  turbulence 
proper  has  set  in  when  each  new  family  of  increasing  un¬ 
stable  modes  is  able,  after  having  reached  a  ct  rtain  finite 
size,  to  cavse  the  generation  c£  a  subsequent  family  of  un¬ 
stable  modes,  etc.  The  first  ones  in  this  chain  of  consecu¬ 
tive  unstable  perturbations  may  be  said  to  belong  to  the  low 
frequency  range,  the  subsequent  ones,  which  generate  each 
other  by  an  inductive  process,  form  the  intermediate  fre¬ 
quency  range. 

"Consider  a  fluid,  which  will,  for  the  sake  of  simplicity, 
be  considered  as  (strictly)  non-viscous  and  (strictly)  iri- 
finitoly  (electrically)  conductive.  If  an  electromagnetic 
field  is  present,  then  the  magnetic  lines  of  force  are,  as 
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is  well  known,  'frozen'  to  the  fluid — that  is,  for  all  mo¬ 
tions  of  the  fluid  each  line  of  force  will  move  strictly  with 
the  material  particles  with  which  it  coincided  originally. 

If  the  fluid  executes  a  complicated  motion,  then  these  lines 
of  force  will  assume  very  involved  shapes,  they  will  there¬ 
fore  produce  large  gradients,  that  is,  large  magnetic  forces, 
and  hence  large  electromagnetic  energies.  If  the  fluid  is 
now  allowed  a  small  viscosity,  and  the  electromagnetic 
field  is  at  first  weak,  then  the  fluid's  motion  will  be  es¬ 
sentially  turbulent,  and  the  electromagnetic  energies  will 
be  greatly  increased.  This  mechanism  will  continue  to 
play,  as  long  as  the  electromagnetic  field  is  small.  Thus  it 
is  anamplifying  mechanism,  which  continues  to  operate 
until  it  has  amplified  itself  out  of  the  linear  (the  amplify¬ 
ing)  range"— that  is,  until  the  electromagnetic  phenomena 
have  risen  to  the  same  order  of  importance  as  the  purely 
hydrodynamical  ones. 

"Let  the  fluid  now  possess  a  large,  but  not  infinite,  con¬ 
ductivity  (i.  e.  ,  a  small,  but  not  zero,  resistivity).  Then 
the  phenomena  will  occur  essentially  as  desc’-ibed  above, 
except  that  the  magnetic  lines  of  force  do  now  possess  a 
small  mobility  with  respect  to  matter.  This  mobility, 
no  matter  how  small,  will  certainly  become  important  in 
the  regions  of  sufficiently  high  involvement  of  the  magnetic 
lines  of  force  —  that  is,  where  they  have  become  sufficiently 
tortuous,  and  sufficiently  closely  pressed  together.  Here 
the  finite  conductivity  (i.  e,  ,  the  non -zero  resistivity)  will 
cause  an  energy  dissipation. 

"Note  the  similarity  between  this  mechanism  and  that  one 
of  dissipation  by  viscosity;  Any  non-zero  viscosity  co¬ 
efficient,  no  matter  how  small,  must  be  able  to  take  care 
of  any  given  turbulent  dissipation  of  energy,  by  failing  to 
inhibit  the  formation  of  sufficiently  small  eddies,  and  then 
working  on  the  correspondingly  high  velocity  gradients 
which  are  present  in  these.  Similarly,  any  finite  conductiv¬ 
ity,  no  matter  how  large  (i.  e,  ,  any  non-zero  resistivity, 
no  matter  how  small),  may  be  able  to  take  care  of  a  given 
dissipation  of  energy,  by  failing  to  inhibit  the  formation 
of  sufficiently  involved  and  sufficiently  closely  compressed 
magnetic  lines  of  force,  and  then  working  on  the  high  fields 
created  by  these  and  by  their  motion.  Thus  there  are 
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'electromagnetic  eddies  very  much  like  the  purely  hydro- 
dynamical  ones,  and  conductivity  (or  rather  its  reciprocal, 
resistivity)  may  be  able  to  bear  to  turbulence  very  near  the 
same  relationship  as  viscosity. 

"From  the  point  of  view  of  theoretical  physics,  turbulence 
is  the  first  clearcut  instance  calling  for  a  new  form  of 
statistical  mechanics:  As  the  earlier  discussions  have 
shown,  it  requires  the  establishing  of  general  statistical 
laws  for  energy  transport  phenomena  in  the  Fourier-trauis- 
form  space,  that  is,  of  the  adjustment  of  the  degrees  of 
freedom  of  a  system,  which  is  endowed  with  many  degrees 
of  freedom,  to  a  flow  of  energy  which  enters  it  in  one  set 
of  frequencies  and  leaves  it  in  another.  (....)  The 
existing  theories  (especially  that  one  of  Kolmogroff 
Onsager  and  Heisenberg)  suffice  to  show  that  these  laws 
will  differ  essentially  from  those  of  classical  (Maxwell- 
Boltzmann-Gibbsian)  statistical  mechanics.  Thus  it  is 
certain  that  the  law  of  equipartition  of  energy  between  all 
degrees  of  freedom,  which  is  valid  in  the  latter,  is  replaced 
by  something  altogether  different  in  the  former.  " 

In  addition  to  the  standard  references  listed  above,  we  have  relied 
upon  Turbulent  Flows  and  Heat  Transfer,  editor;  C.  C.  Lin,  Princeton 
University  Press,  for  the  literature  up  to  1959.  The  pertinent  sections 
of  this  book  are  Chapters  2  and  5  of  Part  B  and  Chapters  1  through 
5  of  Part  C.  For  references  up  to  1962,  we  have  relied  heavily  on 
International  Symposium  on  Fundamental  Problems  in  Turbulence 
and  their  Relation  to  Geophysics,  Journal  of  Geophysical  Research, 
vol  67,  no,  8,  July  1962.  Of  all  the  papers  presented  in  this  symposium 
titles  and  abstracts  of  those  most  pertinent  are  the  following: 

1.  Turbulence  in  the  Presence  of  a  Vertical  Body  Force  and  Temper¬ 
ature  Gradient,  by  Robert  G,  Deissler, 

Two-point  correlation  equations  which  include  the  effects  of  a 
xmiform  temperature  gradient  and  body  force  are  constructed  from  the 
Navier-StokeSv  heat-transfer,  and  continuity  equations.  A  solution 
is  obtained  b-  onverting  the  correlation  equations  to  spectral  form 
and  assumi  g  that  the  turbulence  is  sufficiently  weak  for  triple  cor¬ 
relations  to  be  neglected.  It  is  shown  that,  the  turbulence  decays  with 
time,  although  the  rate  of  decay  is  altered  by  buoyance  effects  caused 
by  the  body  force  and  temperature  gradient.  The  buoyancy  forces  can 
either  exti  .ct  energy  from  the  turbulent  field  or  feed  into  it,  depending 
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on  the  directions  of  the  body  force  and  temperature  gradient.  Spectra 
are  calculated  for  the  turbulent  energy  and  for  the  various  terms  in 
the  turbulent  energy  equation  as  well  as  for  the  temperature  fluctuations 
and  turbulent  heat  transfer.  For  fluids  with  Prandtl  numbers  less  than 
1  the  buoyancy  forces  act  mainly  on  the  large  eddies,  whereas  for 
higher  Prandtl  numbers  they  can  act  on  the  smaller  ones.  When  the 
buoyancy  forces  are  stabilizing,  the  turbulence  can  cause  heat  to  flow 
against  the  temperature  gradient  for  certain  values  of  the  parameters. 
For  making  the  calculations,  it  is  assumed  that  the  turbulence  is 
initially  isotropic  and  the  temperature  fluctuations  initially  zero. 

2.  Some  Mathematical  Models  Generalizing  the  Model  of  Homogeneous 
and  Isotropic  Turbulence,  by  A.  M.  Yaglom. 

Homogeneous  and  isotropic  turbulence  is  an  example  of  a  system  of 
random  fields  invariant  with  respect  to  a  group  of  motions.  Along 
with  homogeneous  and  isotropic  fields,  locally  homogeneous  and  locally 
isotropic  ones  play  an  important  role  in  turbulence  theory;  such  local 
fields  may  also  have  an  accurate  mathematical  definition.  The  random 
fields  invarirnt  with  respect  to  groups  of  transformations  different 
from  a  group  of  Euclidean  motions  can  also  be  considered;  the  'spectral 
representation'  of  such  a  field  and  of  a  corresponding  correlation  func¬ 
tion  often  has  an  unusual  form,  although  its  sense  remains  the  same. 

The  algebraic  theory  of  grovxp  representations  gives  the  general  method 
of  obtaining  the  spectral  representation  for  the  fields.  The  examples 
of  homogeneous  random  fields  on  a  sphere  and  fields  in  a  semiplane 
invariant  with  resnect  to  all  similarity  transformations  present  inter¬ 
esting  examples  of  random  fields  invariant  with  respect  to  'motions' 
of  special  type  which  might  be  of  some  importance  for  turbulence  theory. 

Some  other  related  papers  by  Yaglom  follow: 

Homogeneous  and  Isotropic  Turbulence  in  a  Viscous  Compressible 
Fluid,  Izvest.  Akad.  Nauk  SSSR,  Ser.  Geograph,  and  Geophys.  , 

12(6),  1948;  see  also  German  translation  of  this  paper  in  Statistische 
Theorie  der  Turbulenz,  Akademie-Verlag,  Berlin,  1958, 

Some  Classes  of  Random  Fields  in  n-Dimensional  Space  Related 
to  Stationary  Random  Processes,  Theory  of  Probability  and  Its 
Application,  2(3),  1957. 

Positive  Definite  Functions  and  Homogeneous  Random  Fields  on 
Groups  and  Homogeneous  Spaces,  Doklady  Akad.  Nauk  SSSR, 

135(6),  I960. 
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Second-Order  Homogeneous  Random  Fields,  Proc,  4th  Berkeley 
Symposium  on  Math,  Statistics  and  Probability,  2,  Berkeley- Los 
Angeles,  1961, 

3.  Energy  Transfer  in  an  Isotropic  Turbulent  Flow,  by  Yoshimitsu 
Ogurit. 

This  paper  examines  the  dynamic  consequence  of  the  hypothesis  that 
fourth-order  mean  values  of  the  fluctuating  velocity  components  are 
related  to  second- order  mean  values  as  they  would  be  for  a  normal 
joint-probability  distribution.  The  equations  derived  by  Tatsumi  for 
isotropic  turbulence  on  the  basis  of  this  hypothesis  are  integrated 
numerically  as  an  initial  value  problem  for  an  inviscid  fluid.  The 
most  remarkable  feature  revealed  by  the  computation  is  that  the 
energy  spectrum  function  becomes  negative  during  the  course  of  time 
in  certain  regions  of  wave-number  space.  This  situation  is  similar 
to  the  result  obtained  previously  for  two-dimensional  turbulence. 
Truncation  errors  that  arise  from  finite-difference  approximations  in 
numerical  integration  are  examined.  It  is  tentatively  concluded  that 
this  uuphysical  negative  energy  is  not  generated  by  the  truncation  errors 
but  is  the  consequence  of  the  quasi -normality  hypothesis, 

Reference 

Tatsumi,  T,  ,  The  Theory  of  Decay  Process  of  Incompressible,  Iso¬ 
tropic  Turbulence,  Proc.  Roy.  Soc.  London,  A,  239»  16-45, 

1957. 

4.  Perturbation  Analysis  of  the  Navier-Stokes  Equations  in  Lagrangian 
Form  with  Selected  Linear  Solutions,  by  Willard  J.  Pierson. 

The  Navier-Stokes  equations  for  incompressible  flow  in  their  Lagrangian 
form  are  taken  as  a  starting  point.  A  perturbation  technique  is  then 
used  to  obtain  first-  and  second-order  sets  of  equations,  and  the  gen¬ 
eral  procedure  for  solving  the  equations  to  any  order  is  given.  The 
first-order  equations  yield  interesting  two-  and  three-dimensional 
motions  that  have  some  of  the  properties  of  'stirring,  '  'eddies,  '  and 
'turbulence;'  it  is  suggested  that  various  problems  in  turbulent  motion 
might  possibly  be  re-examined  by  means  of  these  equations. 
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In  addition,  we  have  made  a  systematic  search  of  Mathematical  Reviews, 
Science  Abstracts,  and  the  DDC  abstracts.  The  papers  bearing  on  the 
problem  at  hand  with  abstracts  when  available  are  listed  below. 

Elliott,  L.  A.  ,  Similarity  Methods  in  Radiation  Hydrodynamics,  Pioc. 
Roy,  Soc.  ,  vol  258,  A.  (25  October  I960). 

The  extent  to  which  problems  in  radiation  hydrodynamics  may  be 
tackled  by  similarity  methods  is  investigated  for  the  case  of  spherical 
symmetry.  In  particular  the  problem  of  an  intense  explosion  is 
considered;  by  the  introduction  of  radiative  heat  flux  the  singularity 
in  a  solution  by  Taylor  (1950)  is  removed. 

Chandrasekhar,  S.  ,  On  Turbulence  Caused  by  Thermal  Instability  , 

Proc.  Roy,  Soc.  ,  vol  244,  A,  884  (11  March  1952). 

In  this  paper  a  statistical  theory  of  turbulence  in  incompressible 
fluid  caused  by  the  joint  effects  of  gravity,  and  thermal  instability, 
is  developed.  The  mathematical  theory  is  based  on  the  equations 
of  continuity  and  heat  conduction  and  the  Boussinesq  form  of  the 
equations  of  motion  in  which  the  variations  of  density  (resulting  from 
the  variations  in  temperature)  are  taken  into  account  only  in  so  far 
as  they  modify  the  action  of  gravity.  By  restricting  oneself  to  a 
portion  of  the  fluid  far  from  the  bounding  surfaces  one  can  treat 
the  turbulence  as  approximately  homogeneous  and  axisymmetric 
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and  use  the  theory  of  axisymmetric  vectors  and  tensors  recently 
dev-^eloped  by  the  author  (Chandrasekhar  1950a).  A  number  of  correla¬ 
tions  between  the  various  field  quantities  (such  as  the  velocity  compo¬ 
nents,  fluctuations  in  temperature,  etc,  )  at  two  different  points  in 
the  medium  are  defined;  and  a  closed  system  of  equations  for  the 
defining  scalars  are  derived  for  the  case  when  the  non-linear  terms 
in  the  equations  of  motion  and  beat  conduction  can  be  neglected  and 
a  constant  mean  adverse  temperature  gradient  is  maintained.  Under 
stationary  conditions  when  the  time  derivatives  of  the  various  cor¬ 
relations  are  zero,  there  is  an  exact  balance  between  the  dissipation 
of  kinetic  energy  by  viscosity  and  the  liberation  of  potential  energy 
by  gravity.. 


A  fundamental  set  of  .solutions  of  the  equations  governing  stationary 
turbulence  is  obtained;  these  solutions,  varying  periodically  in  the 
vertical  direction,  enable  a  generalized  Fourier  analysis  of  the 
various  correlation  functions.  According  to  these  solutions,  a  Fourier 
analysis  of  correlations  such  as  ~  ~  T  of  the  vertical  velocities 


at  two  points  directly  above  one  another  and  separated  by  a  distance 
z  ,  cannot  include  wavelengths  less  than  a  certain  minimum  value 
depending  on  ^ne  physical  parameters  and  on  the  temperature  gradi>.nt 
maintained.  We  may  thus  speak  of  a  smallest  size  for  the  eddies. 
Further,  it  appears  that  the  field  of  turbulence  can  be  analyzed  into 
two  modes  characterized  by  the  kinetic  energy  being  confined,  princi¬ 
pally,  to  the  vertical  or  to  the  horizontal  direction. 


Mcore,  D,  W,  ,  The  Boundary  Layer  on  a  Spherical  Gas  Bubble,  Journal 
of  Fluid  Mech,  ,  vol  16, 

The  equations  governing  the  boundary  layer  on  a  spherical  gas  bubble 
rising  steadily  through  liquid  of  small  viscosity  are  derived.  These 
equations  are  linear  and  are  solved  in  closed  form.  The  boundary 
layer  separates  at  the  rear  stagnation  point  of  the  bubble  to  form  a 
thin  wake,  whose  structure  is  determined.  Thus  the  drag  force  car 
be  calculated  from  the  momentum  defect.  The  value  obtained  is 
l2naU^,  where  a  is  the  bubble  radius  and  U  the  terminal  velocity, 
and  this  agrees  with  the  result  of  Levich  (1949)  who  argued  from  the 
viscous  dissipation  in  the  potential  flow  round  the  bubble.  The  next 
term  in  an  expansion  of  the  drag  in  descending  fractional  powers  of  R 
is  found  and  the  results  compared  with  experiment. 

Pao,  Yih-Ho,  Growth  of  a  Weak  Magnetic  Field  in  a  Turbulent  Conduct¬ 
ing  Fluid  with  I.arge  Magnetic  Prandtl  Number,  The  Physics  of 
Fluids,  vol  6,  no.  5. 
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Tiic  behavior  of  a  weak  fluctuating  magnetic  field  in  a  turbulent  elec¬ 
trically  conducting  fluid  with  a  large  magnetic  Prandtl  nun.ber  is  re¬ 
examined,  without  using  vorticity  analogy.  For  a  large  magnetic 
Prandli  number,  the  characteristic  length  of  the  small  eddies  of 
turbulence  is  n.u-  h  larger  than  that  of  the  small  loops  of  the  magnetic 
field.  By  inferring  that  the  small  eddies  of  turbulence  are  mainly 
resocnsible  fot  the  change  of  the  fluctuation  intensity  of  the  magnetic 
field  and  by  utilizing  Batchelor's  mixing  approach,  visualizing  the 
turbulent  velocity  field  as  effectively  a  rather  persistent  unifoi  m  strain¬ 
ing  motion  for  the  small-scale  variation  of  magnetic  field,  it  was 
found  that  the  fluctuation  intensity  of  the  magnetic  field  gro\,s  due  to 
turbulent  motion  for  very  large  magnetic  Prandtl  numbers  ^{v/X)  >  100, 
say  j  and  approaches  a  constant  value  asymptotically.  The  rate  of 
growth  is  found  to  be  a  function  of  v/X  and  T  (the  product  of  time  and  the 
turbulent  straining  rate). 

Chandrasekhar,  S.  ,  The  Onset  of  Convection  by  Thermal  Instability  in 
Spherical  Shells,  Philosophy  Magazine,  ser  7,  vol  44,  no,  350, 
March  1953. 

I  !  this  paper  the  problem  of  the  thermal  instability  of  an  incompressible 
solivie  consisting  of  an  inviscid  core  and  a  viscous  mantle  is  considered, 
and  it  is  sl^own  that  the  pattern  of  convection  which  sets  in,  at  marginal 
stability,  in  the  mantle  shifts  to  harmonics  of  the  higher  orders  as 
the  thickness  of  the  mantle  decreases.  Thus,  when  the  mantle  extends 
to  a  depth  of  half  the  radius  of  the  sphere,  the  harmonics  of  orders 
three  and  four  set  in  about  simultaneously,  while  the  harmonic  of 
order  five  follows  very  soon  afterwards.  The  bearing  of  this  result 
on  the  problem  of  convection  in  the  earth's  mantle  and  of  the  inter¬ 
pretation  of  the  earth's  topographic  features  i?  indicated. 

Chandrasekhar,  S,  ,  The  Decay  of  Axisymmei.ric  Turbulence,  Proc. 

Roy  Soc,  ,  vol  Z03  (10  October  1950). 

In  this  paper  the  decay  of  axisymmetric  turbuler:e  is  investigated. 
Explicit  solutions,  appropriate  for  the  final  period  of  decay,  are 
obtained:  these  solutions  are  in  agreement  with  Batchelor's  gener  ?! 
results  on  the  role  of  the  big  eddies  in  homogeneou.s  turbulence.  It 
is  further  shown  that,  compati'ole  with  the  equations  which  have  been 
derived  for  axisymmetric  turbclence,  a  state  of  tuibulonce  exists  which 
may  be  pictured  as  a  superposition  of  two  non-interacting  fields  which 
are  isotropic  and  axisymmetric,  respectively, 

Tatsumij  T<  ,  The  Theory  of  Decay  Process  of  Inr ompre.ssible  Iso¬ 
tropic  Tu.i'bulence ,  Proc,  Roy  Soc.  239A,  p  16, 


56 


APPENDIX  B 


One  of  tiie  most  serious  difficulties  in  the  theory  of  homogeneous 
turbulence  is  the  indeterminancy  of  the  equation  for  the  velocity  cor¬ 
relation  function  of  any  order,  each  involving  the  correlation  of  higher- 
by-one  order.  In  the  present  paper  this  difficv.i»y  is  resolved  by 
treating  the  two  dynamical  equations  for  the  second-  and  third-order 
velocity  correlations,  and  by  introducing  the  assumption  of  the  zero 
fourth-order  cumulant  of  the  velocity  field  which  yields  a  relationship 
between  the  fourth-  and  second-order  velocity  correlations.  Actual 
calculation,  however,  is  carried  out  in  the  wave-number  space,  and 
a  pai”  of  simultaneous  equations  for  the  energy  spectrum  function 
are  derived  in  Part  I . 

Another  difficulty  of  the  subject  arises  from  the  present  lack,  of 
knowledge  about  the  initial  state  of  turbulence.  In  Part  II,  some 
probable  initial  conditions  for  the  energy  spectrum  are  examined, 
among  which  the  initial  spectrum  of  single-line  type  is  chosen  as  the 
most  suitable  for  the  present  problem  and  its  dynamical  consequences 
are  fully  discussed,  The  power-series  solution  for  the  initial  spectrum 
as  well  as  the  energy  decay  law  due  to  it  are  computed  and  compared 
with  experimental  data.  It  is  found  that  the  solution,  in  so  far  as  the 
approximate  expression  calculated  in  the  present  paper  is  concerned, 
corresponds  to  the  earlier  initial  period  of  decay,  A  solution  which 
would  be  essentially  in  agreement  with  experiments  is  expected  to  be 
given  by  extending  the  present  solution  to  the  further  developed  stage 
of  decay, 

Morton,  B,  R,  ,  Sir  Geoffrey  Taylor,  and  J,  S,  Turner,  Turbulent 
Gravitati  al  Convection  from  Maintained  and  Instantaneous 
Sources,  r  roc,  Roy,  Soc,  ,  vol  234,  A,  (24  January  1956), 

Theories  of  convection  from  maintained  and  instantaneous  sources  of 
buoyancy  are  developed,  using  methods  which  are  applicable  to  strat¬ 
ified  body  fluids  with  any  variation  of  density  with  height;  detailed 
solutions  have  been  presented  for  the  case  of  a  stably  stratified  fluid 
with  a  linear  density  gradient.  The  three  main  assumptions  involved 
are: 


1,  That  the  profiles  of  vertical  velocity  and  buoyancy  are  similar 
at  all  heights 

2,  That  the  rate  of  entrainment  of  fluid  at  any  he:  jht  is  pro¬ 
portional  to  a  characteristic  velocity  at  that  hei  ght 

3,  That  the  fluids  are  incompressible  and  do  not  change 
volume  on  mixing,  and  that  local  variations  in  density  through¬ 
out  the  motion  are  small  compared  to  some  reference  density. 
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The  governing  equations  are  derived  in  rion-dimensiunal  forni  iron, 
the  conditions  of  conservation  of  volume,  moireptum  and  buoyancy,  ami 
a  numerical  solution  is  obtained  for  the  case  of  the  maintained  source. 
This  leads  to  a  prediction  of  the  final  height  to  which  a  plumf*  of  light 
fluid  will  rise  in  a  stably  stratified  fluid.  Estimates  of  the  constant 
governing  the  rate  of  entrainment  are  made  by  comparing  t'ne  tiuuiiy 
with  some  previous  results  in  uniform  fluids,  .rnd  with  tne  results  of 
new  experiments  carried  out  in  a  stratified  salt  solution. 

For  the  case  of  an  instantaneous  source  of  buoyancy  there  is  an  exact 
solution;  the  entrainment  constant  is  again  estimated  from  laboratory 
results  for  a  stratified  fluid. 

Finally,  the  analysis  is  applied  to  the  (compressible)  atmosphere,  by 
making  the  customary  substitution  of  potential  temperature  for  tem¬ 
perature.  Predictions  are  made  of  the  height  to  which  smoke  plumes 
from  typical  sources  of  heat  should  rise  in  a  still,  stably  stratified 
atmosphere  under  various  conditions. 

Chandrasekhar,  S.  ,  Hydromagnetic  Turbulence.  I.  A  Deductive  Theory. 

In  this  paper  .  deductive  theory  of  turbulence  recently  described  by 
the  writer  (Chandrasekhar  1955a)  is  extended  to  hydromagnetics.  By 
consider!  ig  stationary  turbulence  a*'.d  making  statistical  hypotheses 
of  the  same  general  character  as  in  the  hydrodynamic al  theory,  a  pair 
of  differential  equations  are  derived  for  the  scalars  defining  the  iso¬ 
tropic  tensors  describing  the  correlation  in  the  velocities  and  the 
intensities  of  the  magnetic  field  at  two  different  ooints  and  at  two 
different  times. 

Batchelor,  G.  K.  ,  The  Theory  of  Axisymmetric  Turbulence,  Proc. 

Roy,  Soc.  ,  vol  186.  A. 

This  paper  discusses  a  type  of  turbulence  in  a  uniform  stream  whicii 
is  next  to  isotropic  turbulence  in  order  of  simplicity.  Instead  of 
spherical  symmetry,  or  isotropy,  axially  symmetrical  turbulence 
possesses  symmetry  about  an  axis  which  in  practice  is  usually  the 
direction  of  mean  flow.  The  analysis  is  developed  with  the  aid  of 
invariant  theory,  as  suggested  by  a  previous  paper  by  Robertson. 

The  form  of  the  fundamental  velocity  correlation  is  obtained,  and 
scales  of  axisymmetric  tu;.-bulence  are  defined. 

The  results  of  greatest  practical  interest  concern  the  time  rates  of 
change  of  the  mean  squares  of  the  lateral  and  longitudinal  velocity 
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components.  The  rates  of  change  involve  two  terms,  the  fir.st  represent¬ 
ing  vTscous  dissipation,  and  the  second  representing  a  transfer  of 
energy  from  one  component  to  the  other  due  to  the  finite  correlation 
between  the  velocity  and  pressur  e  at  neighbouring  points.  The  effect 
of  the  velocity-pressure  correlation  is  to  bring  the  two  velocity  com¬ 
ponents  towards  equality,  while  the  effect  of  the  viscous  dissipation 
will  only  be  towards  equality  if  an  inequality  between  the  curvatures  at 
the  origin  of  two  particular  velocity  correlation  coefficient  curves, 
both  of  which  are  measurable,  is  obeyed.  The  rates  of  change  of  the 
mean  squares  of  the  vorticity  components  are  also  obtained. 

Chandrasekhar,  S.  ,  The  Thermal  Instability  of  a  Fluid  Sphere  Heated 
Within,  Philosophy  Magazine,  vol  December  1952. 

In  this  paper  the  problem  of  the  therma  instability  of  an  incompressible 
fluid  sphere  heated  within  and  in  equilib  ium  under  its  own  gravitation 
is  considered.  A  general  disturbance  is  analysed  into  modes  in  terms 
of  spherical  harmonics  of  various  orders,  1  ,  and  the  criterion  for  the 
onset  of  convection  for  the  first  fifteen  modes  .s  found  both  when 
the  bounding  surface  is  free  and  when  it  is  rigid;  and  it  is  shown  that 
in  both  cases  the  mode  1  =  1  is  the  first  to  be  excited 

Walters,  J.  K.  and  J.  F.  Davidson,  The  Initial  Motion  of  a  Gas  Bubble 
Formed  in  an  Inviscid  Liquid,  Part  1,  The  Two-Dimensional 
Bubble,  J.  of  Fluid  Mechanics,  vol  12,  no.  408,  March  1952. 

The  paper  deals  with  the  initial  motion  of  a  two-dimensional  bubble 
starting  from  rest  in  the  form  of  a  cylinder  with  its  axis  horizontal. 

The  theory  is  based  on  the  assumptions  of  irrotational  motion  in  the 
liquid  round  the  bubble,  constant  pressure  within  the  bubble,  and 
small  displacements  from  the  cylindrical  form.  This  theory  predicts 
that  the  bubble  should  rise  with  the  acceleration  of  gravity,  over  a 
distance  of  at  least  the  initial  bubble  radius,  and  that  a  tongue  of 
liquid  should  be  projected  up  from  the  base  of  the  bubble  into  its  interior. 
These  predictions  are  confirmed  by  experiments  which  also  show  how 
the  vorticity  necessary  for  steady  motion  in  the  spherical-cap  form 
is  generated  by  the  detachment  of  two  small  bubbles  from  the  back 
of  the  main  bubble. 
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APPENDIX  C 

EVALUATION  OF  INTEGRALS 


The  inversion  in  Equation  23  follows  from 
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The  last  integral  may  be  evaluated  by  integrating  with  respect  to  r 
both  sides  of  the  following  relation 


,,  /tT  -r^/4t 
e  coskrdk  =  -  e 
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this  gives 


r  -k  t  sin  k  r  ,, 
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The  quantity  erf{  - j  is  the  error  function  given  by 
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Assumptions  leading  to  the  ev 
which  appear  in  the  theory  are  examine 
the  role  of  the  principal  convective  age 
Navie r-Stokes  equations  to  linear  form 
moving  frame  of  reference.  The  heat 
simplified  by  the  same  linearization  pr 
tions  in  wave-vector  space  and  subseqi 
initial  conditions  provides  solutions  wh 
actual  flow  situation  found  in  weapons  t 
Qualitative  and  quantitative  features  o 
velocity  distribution,  rise  rate,  and  te 
A  section  on  conclusions  and 
and  results  which  have  come  from  the 
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theory  of  turbulent  flown  (U) 
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.  The  method  is  based  on  the  work  of 
sonents  of  Fourier  spectra  of  velocity  oara- 
vith  compression  waves  (i r rotational  flow) 
he  wave  vector  are  associated  with  eddy 
systems  in  pressure  equilibrium  w'ith  an 
formalism  takes  on  considerable  simplicity 
otion.  (U) 

aluation  of  certain  convolution  integrals 
id  and  show'n  to  involve  the  largest  eddies  in 
nt.  This  assumption  converts  the  nonlinear 
with  solutions  conveniently  described  in  a 
transfer  and  continuity  equations  are  also 
ocess.  Simultaneous  solutions  of  the  equa- 
lent  inversion  to  physical  space  with  certain 
ich  have  considerable  similarity  to  the 
ests  such  as  vortex  ring  (torus)  formation. 

L  the  predicted  motion,  for  example, 
mperature  are  discussed.  (U) 
recommendations  assembles  the  main  ideas 
study  and  suggests  the  direction  of  future 
problem  and  with  respect  to  the  general 
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